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I-        INTRODUCTION 

A.        GENERAL 

The    basis    for  hrlicoptsr   rotor    analysis   was    d9V9lop9d    in 
the    early    1920«s   when   Slauert      expended   propeller    -heory    to 
the    special   case   of    rotating    wings.         Since    that   tiine,      the 
development    of      the    digital     computer   has      permitted    jnany 
improvements   en   Glauert's    analysis. 

In  order    to      develop    a  method   for      predicting   the    total 
rotcr   power   required    in    forward    flight,      it   is   necessary   to 
develop    a      method  that   accurately    predicts      rotor    iynanics. 
The    prediction      of   rotor    dynamics      in    forward    flight      is    a 
complex   one.      Typically,    a   helicopter    rotor    blade    encounters 
a    flow   environment   which    changes      rapidly    as    it    moves    around 
in    azimuth. 

In   forward   flight,    rotor    blade   sections    are    suoiected    to 
azimuthai    variations    in      net    only  angle    of      attac!<:    but    also 
:iach    number.  As   a    result,         comprehensive    perforniance 

analysis    of    a    helicopter    is  much      more    involved    than    that    of 
a    c  o  n  ve  n  t  i  on  a  1   a  ire  r  a  ft . 

Recent  helicopter  design  trends  have  bee^n  in  the 
direction  of  increasing  the  maximum  forward  velocity 
possible    as    well      as    higher  blade   tip      speeds.         Since    the 
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speed  of    the    helicoprsr    is  added    -o    -he    speed    of    ro-=.-iDr.    3f 
rhe    advancing    blade,      the    highest  relati/e    7eloci-i=s    occur 
a-    the  tip    of    the  advancing      blade.         i^hen    the    sec-i::n    Mach 
number  of    the    blade    tip    exceeds      the   critical    Mach   number   of 
the    airfoil,    ccmprsssibility  effects  result.         Thsse   effects 
include    a      large   iicrease    in      profile    drag   and      change    in 
pitching      moments,      therefore      creating   additional      power 
re  quirsments. 

In   all    current   helicopters   there      is    a    tendency    for    the 
retreating    blade    to    stall.      Just    as   the    stall    of    an   airplane 
wing    limits   the    low    speed      perfomance    of    an   airplane,       the 
stall  of      a    rotor      limits    the    high      speed    potential      of    a 
helicopter.         The  relative   velocity    of    the    retreating    blade 
decreases    as      forward   speed      increases.         However,         the 
retreating    blades  aust      produce    the   same    amount      of   lift    as 
the    advancing   blade.       Therefore,    as   the   relative    v<=locity    of 
the    the    retreating   blade    decreases    with   forward    speed,       the 
blade  angle   of      attaclc   must  be   increased      to    equalize    lift 
throughout    the   the   rotor    disc    area.       As   this    angle    increases 
the    blade    will    eventually    stall    at    some    forward   speed.         For 
these  reasons,      stall   and      compressibility    thus    combine    to 
limit  the    perfcrmance   fo    a  helicopter   at    higher    speeds. 
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For    efficient   iasign    of      a   helicoptar,      ths    h=liccp-=!r 
dssignrr    should      have  the      analytical   iddIs      necessary    -c 
predict    the    performance    of  a   helicopter.         This    repor-    gives 
the    user    of   -he    HPf+1-CV    and  -he    IBM    3033    computer    a   means   of 
calculating   the   total   power      required    i.i    forward,      straight 
and    level    flight. 

3.        0EJ2CTIVE 

The    obj^c-ive      was  to    provide  an      intsrac-ive   computer 
program    which    can  oe    used    during    the    initial    dssign   s-age   -o 
estima-e    the   total      power    required    by   the    main      rc-cr   of   a 
helicopter,      as    well    as       maximum    forward    velocity    possible, 
stall  onset   velocity,    and    bes-    endurance    veloci-y. 

Additionally,      -he      desired    accuracy      chosen   a-      -he 
beginning    of      this   effort    was      to   obr:ain    angle      of    at-.ack 
within   a    plus      or   linus    one  half      degrea,       power    es-ima-es 
within  ten    per   cent    of   the     actual    power    required,      and    to 
obtain   the   needed      accuracy  in    as   short   a      running    time    as 
possible. 
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A.        DESCHIPTICN    OF    PROBLEM 

For    a    helicop-c9r    i.i    level      fligh-,      the    maximani    forward 
speed   is    limited    by    the    power    available    as    well    as    ^tall    and 
compressibility    effac-s.         It    is      therefors   aivaniageous    ro 
develop    a      siapl  9-t  o- use    ccniputer      program    -h  =  -      estimates 
blade     s-all   and     compressibility      effects    on      helicop-er 
performance. 

3.        METHOD    CF    SOLUTION 

In  forward   flight,      tihe      aerodynamic   environment   of   the 
rotor   blade   varies  as   the      rotor    blade    rotates    wi-h   respect 
zo    the    direction   of    flighi.      The    method  chosen    to    obtain    'ne 
objective    accuracy      ignores  any      variable    no-      immediately 
impacting    on    the   accuracy    of     the   derermina-ion    of   angle   of 
ai-ack  and   power   requirements. 

The      me-hcd   utilized      herein    uses      a    combination      of 
momen-*:um      and    blade       =lemen-      theory      to    perform      rotor 
performance   calcula-icns.         This    theory   is    ini-ialiy    used   -o 
datermina    "^ha    inducad,    profile,     and    parasite    power   requirad. 
The   solution    me-hod   then  ut:iliz6s   a    blada    elaaen-    analysis 
-0    predict,      the   cyclic,       collective,      and      inflow    anglas 
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associat'^d    wi-h   nhe    rotor.        These      angiss    ara    -hsr.   a3rd    ir. 
•ths    calculation   of   the   rotor    blade's      angls   of   at-.ack.    a-    -.he 
azimuthal    positions    of   90    and    27  0   dsgreas.         Compr ^ssibilit y 
and    s-all    power   ara    then    estimated    as   a   function   of  angle   of 
attack   and    forward  velocity. 
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III.   D5S3RI  PTION  OF  THE  PR0BL2;i  SOLUTION 

A.   GENERAL 

The   calculat icn/pradiction   of      halicop-ar   perfornince    is 
primarily    a   3tatt9r      of  detsrmiriing   the   power      reqairsd   and 
power  available   over    the    desired    flight   regime.         The    power 
inforraaricn   nay      then    be    used      to   predior.      the    operational 
capabilities   of   the    aircraft. 

For   the      ideal  helicopter      (no    lossesi ,         in    forward, 
straight    and    level   flight,      power  required    can    be    subdivided 
into    five    parts:        induced   power,    P    ,      raquired   to    produce 
rotor  thrust;       profile   power,    ?q  ,      required   to    overcome    -he 
s^in    friction    and  pressure  drag    of    the    .iiain    rotor;      parasite 
power,      ?p,      required   to    overcome  the    parasite   drag  of    the 
helicopter;      compressibility    power,    P,,,    required   to   overcome 
the    increase      in    profile      drag   when      the    tip      :iach   number 
exceeds    the    drag    divergence   Jlach    number   of    the    airfoil;       and 
stall   power,       Pg,      required  to    account    for    the    increase    in 
rotor  torque   as   a  result      of   retreating   blade   stall.         The 
HP41-CV    and      IBM    3033      computer    programs      contained   herein 
provides    a    simple,         quick  means    of   evaluating      the   power 
required    by    an   isolated    rotor. 
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3.        ASSUMPTIONS 

The    niajcr   tasK  in    helicop-er    psrf or.ianca    analysis    is    -hr 
determination   cf    rotor   forces,      angles,      and    power.         Th,= 
nie-hod  chosen    ignores      any  variable   thai    does      no-    iirec-ly 
iapac-   on    the   desired   accuracy.  Therefore,      in    orier    to 

simplify    -he      compit a-icn al  process    used,         -he    followina 
assumptions   are    aade: 

1.  Steady    flow    through  the    rotor    system. 

2.  Small      angle      approximations      are      a      valid 
rspresen- ation   of  the   real  world    phenomena. 

3.  All    blades    considered   are   rectangular    (ncn-ta pered) 
with    only    uniform  twist    possible. 

4.  Hinge  off-set   is    zero       (i.e.,      the   thrust    vector 
passes   through   the   C  .  G. )  . 

5.  The  stall  angle  en  most  helicopter  blades  can  be 
approximated  by  the  angle  that  occurs  at  CLMAX  for  3  2-D 
airfoil. 

6.  The    helicopter    is   trimmed.         This    implies    that    ths 
sun    of   all   the    moments  about    the    center   of    gravity    (C.G.)     is 
zero,      all    forces      are   in    balance,      and      that    no   lateral 
flapping    is    present. 
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7.  At      stall   onsat,         the    valu9      Df    3=c-icn      ir^q 
coefficient    jumps  ip proximatsly    0.08. 

8.  The   rotor  ro-ates    countarclockvise . 

C.  NCMENCLATUHE 

"Standard"    nomenclature  is   used.      Appendix    A   contains    an 
alphabetical    list      defining  all    the    symbols      and   inefaonics 
used    in    the    development    of  this    report.      Appendix    3   coniains 
an    alphabetical   list      of    all    HP4 1-cv   displays      used   in    -he 
HP41-CV    program. 

D.  INITIALIZATION 

It   is    assumed  that   the     initial    design    of   the    helicctsrer 
has    been    completed  and   that      the    helicopters    weigh-    as    well 
as    the   chord,    radias,    tip    velocity,      twist,      zerc-lif-    drag 
coefficient,      and  number    of   blades      of    tae    roior    are    known. 
Finally,       it      is    assumed    that      an   initial      es-imat-ion    of 
equivalent    flat      plate  area   is      available,       and      that    the 
forward    velocity    of    the    helicopter    is    known. 

E.  ROTOR     DYNAf^ICS 

There  are  four  rotor  parameters  whio.i  will  help  expedite 
later  calculations.  The  firs-  of  these  is  the  rotor  advance 
ratio,    y    ;      y    is      the   ratio   of   the      helicopter's    forward 
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v^locizy    ro    the    rotational  Vrlocity.         Tiia    advar.c^ 
be    rsnresentsd   as 


(1) 


y    =      V^:,  Co  s    a 


(dimensiorilsss) 


Applying      small      angle   approximations      tae    advance      ratio 
becomes 


(2) 


U     =       V. 


T 
where , 

V^    =     (    i7    R) 

The      second    dime nsionless      ratio      -hat      needs   to      be 
calculated    is    the  inflow    ratio,    X    .      The    inflow    ra-io    is    the 
ratio  of    the    net    velocity    up   through   the    rotor    systam    to    the 
tip    speed.         For    the    near    hover    case,    \i      <    J .  1 ,    the   Tan  a=   0 
and    the    inflow   ratio    can    be   approximated    by 


(3) 


A 


/CT/2 


In    fcrward      flight,       the      calculation    of      X       requires    the 
determination    cf      the    angle  of      attack    of.    the      rotor    disc, 
Lrttincj    the   anale   of    attack   between      -he    disc    plane   and   th; 
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incoming    free-stream    velocity    be  a-.  ,      and    assuaing    angles    -.z 
be    small,    th.e    inflow    velocity   can   be   caloulazed    using 


(3a) 


X      = 


-  c, 


+    u   ran    OL. 


/A    +1J 


where 

a^      =  -Tan     (Dp/L   )     =    -Tan    (Dp/W) 
Dp    =    (Pp*    550)/   Vp 

The      las-:    dimensionless      parameter    ::ha-      needs   to      be 
calculated   is    solidity.         Solidity   is      the    fraction   of    the 
rotor  disc   area   that    is    composed   of    blades.         For    a  blade   of 
constant    chord    (i.e.  ,    non-tapered)     solidity   can    be   expressed 
as 


{^) 


a     =    (b    *    C)/(?I    *    R) 


Finally,    the    tip    loss    factor,    S^nyr      njust    be    considered. 
The    tip    less    factor    is   used  to      account   for    the    loss    of    lift 
that    a   rotor    blade  experiances   due      tc    flow   from    the    rotor's 
lower  surface    to    its    upper  surface.         The    tip   loss    factor   of 
a    rotor    can    be   approximated  by 


(5) 


\l=    '-' 


-     '/2    *    CT      /    b 
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F.        VELOCITY    CALCULATIONS 

Thars    ara    four  velocities    that      are   of   interes-.         The 
first   of    these   is  the      downwash    velocity    ,      w.         Assaaiinq 
steady  flow  through    the    rotor,      the   dowrii/ash    velocity    can   be 
apprcxiaated    by 


(6) 


w       = 


2P    Aq   Vp 


(f-/s) 


NOTE:      This  equation  is   not    valid   for    small 
values    of    forward    velocity. 


The    second   velocity    that   needs      -o    be    calculated    is    the 
3-all  onset    velocity,      7^  .      The    stall   onset,    velocity    is    the 
velocity    at    which  the   retreating    blade    tip    first    exceeds    the 
static   stall    angle.         Tue    forward  speed   at    which   stall    onset 
is    first    noted   can  be   approximated      by    the    velocity   for    best 
range,    V3R.         This   is    due   to   the    marked   increase    in   profile 
pow9r   required      at  speeds    higher      than    valocity      for    best 
range,         A    typical   set   of    power    curves    for    a      helicopter    are 
shown   in    Figure    3.1. 

The    forward    velocity    for    minimum   ?/7     (i.e.,       best    range 
velocity)       is   easily    found  graphically   on    the   pcwar   required 
curve   as    the    point   where    a      straight    line    tlirough    the    origin 
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POWER  HEQyiHED  VS  YELOCaTY 


0.0        2S.0         SO.O         7S.0        100.0      12S.0      1S0.0      178.0 

FORWARD  VELOCITY  (KT) 


Figure   3.1.       Typical   Helicopter   Power   Curves 

is    tangen-    ro    ^he  curve     (see    Figure    3.2),         Since    --he    power 
curves  are      not    initally    known,       an      analytical    expression 
needs   to    be    developed    which      will   estimate    the    velocity    for 
best    range.         The   follDwing  is      the    derivation    of    ens    such 
ex  or ession . 


21 


As   shown    below,       if    the      parasi-ie    powar   raquirB'!    is    s^- 
squal  rD    the    profile    power   riqairsd   and   tha    squali-y    is    -hen 
solved   for      an   equation    in     which   forward    velocity      is    the 
variable,      the      resul-    is    a      cubic    aquation   in      which    the 
largest    root    defines    the    point      where   tha    profile    power   and 
parasita    power   are   equal. 


P  =  D 

-  p  -0 


1/2PfV3=      1/8  a    CdopA^(1    +    4.25      u^    ) 


(7)        ^3-     (U.25   /Uf)Cdo   Aj^7^    ^2      -     (a/^f)Cdo    k^V 


3  =      0 


Next,      the      largest    root   of      equation    7    needs      to    be 
datermined.         This   can   be    siaply      done   using    the    cubic    roc- 


•quation.       Letting 


p    =  -  (U  .25  a/4f)  Cdo    A    V, 


D   T 


r    =   -(  a  /yf)Cdo    A^V^3 


Sauation    7    can   now   be    written   as 


(7a) 


F  ■    F 


Substituting 
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(    X    -    D/3    ) 

for    Vp  in    equa-ion   7a   yialds   an    equation    of  the    fora 
(7b)  X3    +    aX    -^   b    =   0 

a    =    -1/3   ?2 

b    =    1/27    (2p2    +    27r) 

The    largest   rco-    of    equation   7b    is   then   given    by 


(8) 


V      =       A    +    3 
w  h  er  9 , 


A      = 


+    -v/jli      +  ^ 


V 
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-r 


-  V^ 


nl/3 


-,  1/3 


17 


Bquaticn    8    C2.n.   now   be  used     as   the    initial   approximation    for 
stall   onset    velocity. 

The      third   velocity      which   is      of      interest    is      the 
helicopter's    aaximu.-n    forward    velocity,      7i!lAX.         As    shown    in 


z:i 


Figure    3.2,      the    maximam    forward      v9loci":7    is   aivsn    by    -h' 
in-t  =  rsect ion   of      ths    powar     raguirai  and      pow-sr    a'/ail=i3ls 
carves   for    a    given  gross    weight    and    alricud?    [Hef.    11. 


POWER  REQUIRED  VS  VELOCITY 


MAX  POWER  AVAILABLE 


VMAX 


0.0        25.0         SO.O         7S.0        100.0      128.0      1S0.0 

FORWARD  VELOCITY  (KT) 


Figure    3,2.      VMAX,73E,    and  7S   from    Power    Available   Curves. 

Therfore,      whenever    VF   >    VXAX ,      th€ra    -s    insuff icien *    power 
available    to    sustain    level   flight. 

The    power-liaited    maximum   speed   ."nay    be    estimated    by 
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1/3 
(9)  V  =  "^       (P  -p-D\ 

MAX  '^      AVAIL      "I         ^o'  {ft/s%c) 


Equation    9      can    be   simplified      by   assuming    ihat      tha    pcwer 
reguirad    at   maximum    spsed    is  about    the    same   as    -ha-   required 
at    hover.      Therefore,    assuming 


(10)  (?  -?-?)     =     fP  -o 

^      '  ^      AVAIL        I  o'  ^       HOTTER     "  o 


and    that 


(    P  -    ?   )     =    (   P   )  =  W 

AVAIL  0  I  HOVER  /    2P    A 


equation    9    can   then    be   written   as 

1/3 

where , 


(10b)  v_    =       -/      w 

1  V2    p    A^ 


The  las-,  velocity   that  needs  to  be   calculated  is  the 
best  endurance  velDcity,  VBE.   In  the  normal  operating  range 
thi  total  helicopter  power  can  be  represantsd  by 


(11)      ^T  ='l  *  ^    *  ^ 


Assuming    tha-    the      variation   of    profila    pow^r      wi-h    forward 
velocity    is      nagligibls,       -hs    velccizy    for      besr    sniurancs 
(also   the    best   rate    of  climb    velocity)       ran    be    found.         If 
equation         is      differentiated   with      respect      to      forward 
velocity,    VF,    and'  is    set    equal  to   zero,    it    can    be    seen    that 


(11a) 


or 


(11b) 


W2 


^  I      -^  ^p 


=       3  D  f    7j? 


2    P    A^^BE 


Solving   equation        for  the   best    endurance      velocity,      VB^, 
results    in 


(11c) 


VBE   = 


P     A, 


3f 


(ft/s) 


G.        INITIAL    POWER   CALCULATIONS 

As   forward      velocity    increases,         the    induced      power 
decreases,      the    profile    power    increases    slightly,      and    the 
parasite    pcwer   increases    until    it      becoaes    the    dominant    loss 
at    high    speeds   [ Ref .    2].         For   forward,      straight    and    level 
flight,    the    induced    power    can    be    calculated   by 
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(12) 


h 


2^^ 


7v2 


2  +       1 


Y 


v-j.    /    5  50 


(h?) 


If    rip      losses  are  takan      into   affect    (sjaa-ion      5),      the 
induced    power    r.ow  becciaes 


(12a) 


^I  =    <    ^^\l^    *   "l 
(TL) 


Sir.cs  the    induced      power    required   in   ground      effect   is    lass 
than    that    required   out    of    ground    effect,       equation    12    should 
be    written   as 


(12b) 


(TL+GE) 


wnere. 


GE    =     (-0.1276    *     (h/D)  *    +    0.708    *     (h/D)  3 
-1.4569    *     (h/D)  2    +    1.3432    *     (h/D) 
+    0.5  147) 


The    profile   power  required  is    given   by 


(13)  Pq=   aCdo  p  Aplf^     (1    +    4.25  U2> 


4400 


(h?) 


Finally,    parasite  power    given    by 
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(14) 


^p 


of     VL^ 


¥ 


1100 


(h?) 


H.        ANGLE    CF    ATTACK 

Th8    determination    Df      angla    of    attacic    at      the    aziiiuthai 
positions    cf      90    and      270    degrsss      are    inportant      in    th= 
d  =  terminat  ion    cf    co  oa  press  ibilii/    and   stall    effects.         It    is 
because    cf   tiiis    effec":   that  the    dynamics    of    the    blade    T.otion 
are    important    in   analyzing   helicopter    p3rf ormance. 

The    angle      of  attack      of   the    rotor      is    a      function    of 
radius,      r,      and   azimuthal     position,     ^     .         Figure   3.3 
illustrates    the      sign    convention    used    in      detersiining    blade 
angle   of    attack.      The   angle   of   attack    cf    a    rotor   can    thus    be 
estimated    by. 


(15)  a    {T,ib) 


(15a) 
(15b) 


0  -    Vp3cos  a  +    -rf    -    7  a^ 


fiR   +    VpSin   ^ 


(rads) 


w  h  er  a , 
0 

6 


0„+     9^+  ©Cos  'F  +   esin^  +    K^  6 

U  1  1  2  o 

3q-    a1    Cos  ip  -    b1    Sin ']; 
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Figure    3.3.      Rotor    Angles    in   Longitudinal    Plans. 

To    ietarmine      angla    of   attack,         i-    is   -hsn      necessary    to 
ae-isrmina    the      longitudinal  coliec-ive    a^d    oyclic      =nal9s    0 
and    9      .         This      cin    be      accomDlishad    by       expressing    -he 
coefficient    of   -hrust    as 
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(16) 


wher^ , 


aa     aT1     +     (0  +     K„3)T2    +0^13    +     (Q    -    KD^)7^ 
"^  0  B  T  2  B 


T1     =    .5  (32    ♦    . 5  y  2)  T3    = 

TL 


.25B2<B2     ^     y2) 
TL   TL 


'2    =     (.5a3    +    .5u23    )  T4    = 

TL  TL 


. 5  u   (32    +    .25  y  2  ) 
TL 


Additionally,      the      longitudinal    flapping      coefficient    a1, 
nseds     to    be      de-ermined.  Tha    longit,aiinai      flapping 

coefficient   can    be  written  as 


(17)  a1    =AA11     *     (9+    5eB)Al2    +0A13    *     (0-    K    b1)Al4 

0       B   0  T  2  B 


where , 
D1 


(32 

TL 


-    .5    \i2^ 


All     =    .25  (.5      3^2^  -    y3    /3) 

3    2     D1 
TL 


A 1 3     =        2  y   B     2 

^     TL 


D1 


A12    =         8  y   3, 


TL 


3    D1 


A14     =       3    2+1.  5y2 

TL 
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Assuming    tha-      there    is    no      lateral    flapping    and      -.hat    the 
effect   is    zero    (i.  =  ./    Kg    =  0.0),     aqaations    16    and    17    beccm= 
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(18) 


2    C-Y  =        XT1+0r2+0T3+OTU 
"a  0  T  2 


(19) 


a1 


A   All    +0   A  12    +    0mA  13    +0  AlU 

0  -■-  2 


NOTE: 


Since    tha   analysis    is    cni'/   i< 

th9   azimuthal    Dcsitions    or    90    and''270 


.ocKir.g    a' 


ion 


degrees,   i-    can    be    sean    iroi    eguati. 
15  5  tha*   the    ccn-rib'icion    of    tae 
longitudinal    flapping    coefficient,    a1, 
wilx    always   be  zero   at    these    positions 


Equations  18  and  19  new  represent  a  se-  of  simal-anecus 
equations  in  which  the  only  unknowns  are  rhe  collective  and 
cyclic  angles   and     can   thus   be  dexermined.  Knowing    the 

values  of  the  cyclic  and  collective  angias,  zhe  blade  tip 
angle  of  atiack  at  the  90  and  270  degree  posi-ions  can  be 
estimated    by 


(20) 


^90    =    ^    ^    ®T    ^    ^    ^       ' 


1+y 


(21) 


a 


270 
NOTE: 


0       -    0       +0       + 
0  2  T 


X 
1+   U 


The  angle  of  attack,  a  ,  is  defined 
oositively  if  the  disk  olane  is  nos 
iiD,    see  Figure    3.2. 
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I.        STALL    POWHB 

In  the   helicopTior      stall  normally    star  is    a-      -hs    -ip   of 
the    rstrsating   blade,    sines  the    highest   angles   of   attack    sre 
usually    at   the   blade. tip.        As   the    forward    spead   increases, 
the    stalled   area    of    t h9    rotor    blade   spr9ads    inboard. 

At  the    higher   values    of  y,    the    effects    of    stall   on    power 
raguirsd   are      gr^at    and    therefore      ne^d   to      be    esti.iiat  =d. 
Assuming    a    jump   of   0.08    in  the   valua   of   Cdo      at    stall    onset, 
taat    tha    rotor   area    within  which    tha    blade    stall    exists    is   a 
segment    of    minimum      d imensionles s   radius    X^      and      that    the 
stall   area    is    symmetric    about    'p    =   270    degrees.      Castles    and 
New    found    that   the   effects   of    tip   stall    on    power    required    at 
the    higher   values  of   y    are  large    and   can    be    approxi.-aated    for 
high    speed    flight   by 


<22)         cps  =   i;nr  <'  ->^ ''  ('  -  ''s'    ''^"^ 


where   1^    is   the   nondimens icnal      radius    outboard   of   whici   the 
retreating    blade    is      stalled   [Ref.    3].         The    dimensionless 
radius,    X^    /      can  be    estimated   by   equating    the   section    angle 
of    attack,    at    'b  =  270   degrees,       to    A^IAX   ^  Ref .    a].         Setting 
equation    15    at    'p  =   270   degrees      equal    to    A.1AX    results    in    the 
quadratic    listed    below. 


32 


(23)  (X    -y)    =     0^(X    -y)*  0^  X2- e^y-e    (X   -y)+   X  ^1    =1 

Applying    the   quadratic      formula    to    equation    23      yields    -he 
roots 


(23a)  Xg     =     -Bg    +      /B|   -   4  0,pC< 


20^ 


(23b)  X        =      -B       -      /B2  -   a0    r 

0  J o 1     b 

2  0^ 
wh^re, 

'"='  r     =     "max  -  \*   \ 


(23b)  Cg    =      yr  +  A 


(23s)  B     =  -yQ^-r 


Equation    22    is      satisfacrory    for   most    cases.  Iz    is 

possible    however,      for  -^he   blade    sec-ion    angle   of    a-^tack   -zo 
be    higher    inboard  than   at    the    -ip   crea-iag    =    si-uation    which 
is    ussually   refered    to  as    inboard   stalling    [  Sef .    5], 

For  the  special  case  of  inboard  stalling,  see  Figure 
3.4,  the  incremen-al  stall  power  coefficienr  defined  by 
equation      22    is      too    large     and    needs      to    be      corrected. 


INBOARD  STALLING 


VF 


V'  =  71 


i/  =  3Tr/2/ 


V  =  0 


i>  =  n/Z 


Figura    3.4.       Inboard   Stall   Patsrn. 

Assuming    -he    stalled      region    is    dianiond   shaped      as    shown    in 
Figure   3.5    and   that    the    stalled      area    is    syniine-rio    abcu-    '^    = 

3TT/2,      it       can   be  seen   that     as    X      aopr^aohes    X      ,      the 

0      " '  S 

correction    to   the  incremental    stall    power    cDefficient,       '^pa* 
must    vanish. 


3n 


STALL  REGIONS 


(j,l  >  fx2  >  IJ.3 


VF 


^  =  3tv/2 


V  =  0 


ii  =  tt/2 


Figurs  3.5.      Approximated   Stall    Area. 

Similiarly,      as    -the    average  of      X      and    1^    aoproach^s    uni-y, 

0  S 

-he    value    cf    Cpg   gDes   zo    "he   value   defined    by    equation    22. 

Therefore,      in  order   to  correct      for    tha    oossibiii-^.  7   cf 
ip.  board   stalling,    -he   correction    factor,    !c    ,    is    defined    such 
that 
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(24) 


k      =       1.0 
S 


fD-      -3 


20. 


A     >     1.0 


(25) 


k      =      B    /20    +   X^ 
S         T         S 

1  -  % 


ZQ 


<  1.0 


The    corrected   aquation   for  stall    power   raqiirsd   "hus    becomes 


(26) 


PS  S  PS 


J.        COMPRESSIBILITY    POWER 

Ths    individual   sffects  of      stall   and    coaiprsssibili'-.y    on 
ro-or   profile    power    are    substantial      at    liigh    aivanca    ratios. 
When    both    effects   are   present    the      losses    due  to    each    source 
are    difficult    zo    distinguish.         Therefore,    as    a    helicopter's 
forward   velocity   and      tip   speed    increase,      the      need    for   a 
simple    estimate      of    how      the   compressibility      of   the      air 
influences   the   rotor    performance    is   necoessary. 

In   forward   flight,         the    Mach   number    of      the    advancing 
blade   is    given    by 
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(27)  :i       =       M^^p  (X    ■»■  y  Sin  .j;  ) 


(27a)  M^^p    =    fiR         =7^ 


(27b)  X      =      r/R      =      nondiasasional    radias 


Since  'he  highest  Mach  number  occurs  \z  -h9  lip  of  -h= 
advancing  blade  at  i|j  =90  degrees,  ejaation  27  ::a:i  be 
written    as 

(28)  M        =       M       (1    *  y  ) 

Reference    3.5    showed    -ha*    the      critical   ."lach    number   for    drag 
divergence   can   be  estimated  by 


(29)  M^      =       .^CRO    -  ml    *    a    *  a 

CRIT  90 


Gesssow      and      Crim,         in      their      investigation      of 
compressibility   found      that  the    compressibiliiy      effec*    on 
roior   performance   was   a    rapid      increase    in    -he    profile    pcwe: 
when    the    tip    Mach  number      exceeded   -che    critical    ilach    number 
for    drag    divergence    [Ref.    6].       The    incrsaee    in    profilr    pcwe: 
coefficient    due   to   !*ach    effects    can    be    estimated    bv 
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(30)  Cp(.  =  0.01   AMj3  +    O.K   a:1d>  ' 

where , 

(31)  ^%'    "^       -    "CMT-     ^-"^ 
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IV.       CDNCLUSICNS    AND    RECOMEfi D  ATION  S 

Ths    objactive   of    -his    project   was      -o    proviaa    ar.    ^asy   -.0 
use,      m-sr-active     ccmputer   program   whi^h    could      br   used 
during  the    initial  design    phas?      to   estimaza   the   -oral    power 
required    by  -he    main    rotor     of   a    helicopter.         The   compu-er 
programs    con-ained      herein  provide   results    v/hich      are    well 
wi-hin  the   objective    accuracies     (see    Figures    C.I,    C.2,      C.3) 
and    provide   acceptable      resul-s    as   a   first      cut    estimate    of 
compressibility   and    stall    power    requirements. 

In  the      development    of  the   computational      model,      many 
simplifying   assumptions      were    made    to      ease    the      airount    of 
computation    required.         The  assumption    which    most    impacts    on 
the    accuracy    of    the    program  is      that   of    steady   flow   *-.hrough 
the    rotor. 

The    flow      environment   encountered    by    the      rotor   changes 
rapidly   due   to   the  rate    of     change    of    blade   angle    cf    attack. 
Additionally,      rotor    operation   at      high    idvince    ratios    also 
produces    considerable   radial      flow   along    the    blade      span    . 
The    steady    flow    assumption     ignores    the   time-variant    aspect 
of    rotor    aerodynamics.         The   dynamic      nature    of    tn ?   rotor, 
especially    when      operating   at    Dr      near    the      stall    regime, 
requires    the    application    of      unsteady    aerodynamics,      and    a 
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clo£€   examination  of    how    ths      pitching    oiomants    gsr.  era-:  ri    by 
"he  retreating  blade  stalling  affects 

controiability  [Hef.    7].         it    is      -herefore   recomended    that 
additional      investigations     consider      how      -he      uns-eady 
aerodynamics,    and  pitching  moments   generated,      influence    tn- 
performance   of   the   helicopter. 
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APPENDIX    A 


NCMENCLATURS 


T§r2      Hne mcnic 
a  a 

A  AD 


a-.  Alphas 


^'max        ^''^^ 


a  A9  0 

90 


a  A270 

270 


a1  a1 

All  All 

A  1  2  A  1  2 


Definition 
Slope    of  airfoil 
section    lift   curva. 
Rotor    disc   area. 
Disk   plane   angle   of 
attack. 

The   steady    flow   stall 
angle    of  the   airfoil 
(given    by      a   *    CLMAX   ) 
Angle    of  az-ack   of 
the   advancing    blade 
at  ip  =    90   degrees . 
Angle    of  attack   cf 
the   retreating    blade 
at  ^li    =    270   deg  rees. 
Longitudinal   flapping 
Terni   in    definition   of 
rHSTA2. 
Term   in   definition   of 


Unit; 
rad 


•ad 


raa 


rad 


rad 


rad 


dime  r.sionl  ess 


diiaensionless 


rHETA2. 
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A  1  3  A  1  3 


A  1  U  A  1  U 


Tarm   in    definiton   cf 

rHETA2. 

Tarni   in    definition   of 


iims'isi  onl  r  ss 


din^r.-^ionl^ss 


rHETA2. 


Number    of  blades   on 


iim9r.3ionl<rS3 


8  3ETA0 

0 


Sg  B3 


Btl        btl 


the   main  rotor. 

Main   rotor  coning 

angle. 

Term   in   definition  o. 


^s  • 


Tip   loss  factor. 


rad 


dia^nsionlass 


dimensionleso 


Cdc  CDO 


MAX 


C  C^  C 

PC 


c  c^  s 


Cg  CS 


Mean   chord  of    main 
rotor   blade. 
Main   rotor    coefficieni 
of    drag   at   zero  lift. 
Maxiaum    coeficien-    of 
lift    (2-D). 
Correction   to    powar 
coefficient   due   to 
compressibility   effects. 
Correction   -o    power 
coefficient    due   -o   s-all 
Per  a   in    definition    of    i^ 
'42 


a:..Tiensioniess 


dimensionless 


(iimensionxsss 


dimensionless 


21397-31  Dr.^~5S 


Dp 
D 

All, 


G.  E. 
h 

K 


D 
DP 

DMD 


6  DEL3 

3 


?'?  k 


GAMMA 

GE 
H 

HP 
K3 


Coefficient    of    thrust 

of    the    main    rotcr. 

Rotor    disc   diame-er. 

Parasite  drag    of    the 

helicopter. 

Tern  in    the   definition 

of    compressibility 

p  0  w  er . 

Sate  of   change    of   bladf 

pitch    with   respec-    to 

blade    flapping. 

Squi7al€nt    flat   plate 

area   of    the    helicopter 

in    forward    flight. 

lerm  in    definition   of 

B      and    C     . 

Ground    effect    ratio. 

Height    of  main    rotor 

above    the  ground. 

Horsepower. 

effect . 


diiuensionl-TSS 


dimensionless 


diaeneionless 


dimensionless 


no 


iiaensionless 
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^  LAMB  Ratio    of  the   ne"c  iinier.'^icnl^ss 

velocity   ap    through 

the   rotor  system    to 

the   tip    spe^d. 
!!  MACH  iHach   number    of    rcior  dims-nsiorLl  rss 

blade. 
'"^CRIT      '^C^^T  Critical   Macn    number  iimer.sionless 

of    advancing   blade 

at    ip  =    90  degrees. 

H^  MT  rio    .^ach  number  of  dimensionless 

T 

rotor    balde. 

M^„„^      MCRTO  Cri-ical   Mach    number  dimensicnisss 

CRTO 

for   CI    =  0. 0    . 

y  MU  Sain   rotor    advance  dimensionless 

ratio. 
ml  Ml  Constant   in    definition  ft-lb 

of    critical    Ma ch    number 

(ml    =    0.  113)  . 

Q,  OMEGA  Rotational   veloci-y.  rad/s 

?  ?!l  Power    required    due   to  hp 

C 

compressibility   eff=c-3, 
?  ?I  Induced    power.  hp 


0 


?s 


T1 


T2 


PO 
?? 
PS 


R 

R 

P 

RHO 

a 

SD 

e 

THETA 

THETAO 


theta; 


T1 


T2 


r3    13 


T4 


Profila  power. 

Parasite  Pow=r. 

Power    required    due   to 

stall    effec-^.s. 

Main  rotor  radius. 

Air  density. 

Main   rotor   solidity. 

Ratio    of  ambient    temp- 

cra-ure    to    standard   S9a 

level    -cemperat  ure . 

Main   rotor   collective 

pitch. 

Main   rotor    lateral 

pitch. 

Main  rotor   longitudinal 

cyclic    pitch. 

Term   in    definition   of 

THETAO. 

Term   in    definition   of 

THETAO. 

Term   in   definition   of 

THETAO. 

Term   in    definition   of 


hp 


no 


ho 


slua/ft 


di-iiiensioniess 


.im^nsionless 


rad 


rad 


rad 


dimensionless 


diaensior.l  ess 


i i ma ns i 3 nl ess 


dime nsi unless 
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V  J  VI 

7^  VF 

7  7KT 

VMAX  VJIAX 


sv 


SVEL 

DW 

XS 
XO 


r  H  E  TA  0  . 

Induced    velocity. 

\ircraft   forward   speed. 

Aircraft   forward   speed, 

Saximum    forward  velocity 

possible. 

S-all    onset    veloci-y 

(velocity  at    which    A270 

equal    A  MAX) 

Rotor    tip  speed. 

Speed    of  sound  . 

aotor    dcwnwash    velocity. 

Aircraft   gross    weight. 

Radius    outboard   of    which 

the    main   rotor    is    s^all^d. 

Radius    inboari    of    which 

rotor    blade    stall    may 

be    present    due    to    inflow 

ratio    and   blade    twist. 


ft/s 

f-/s 

kt 


ft/3 
fz/3 
ft/3 

lb 
dimensionless 


dimension less 


U6 


APPENDIX    B 
Hi   llz^l    PROGRAM     DOCU^JNIATION 

A.        GENERAL 

The  HP  41-cv  program  uses  45  program  storage  r=gis-=rs. 
Prior  to  program  initialization,  -h5  handheld  calculaicr 
should   09    sized   "o  46. 

Tha  ccmputsr  program  consisis  of  a  main  program  which 
calculates  the  necessary  dynamic  param^-ars  of  the  main 
rotor  and  several  subroutines  which  calculate  power  required 
as  well  as  the  maximum  forward  velocity,  s-all  onss- 
veloci-y,    and   best   endurance   velocity. 

The  documentation  for  the  HP  4  1-CV  program  is  aevided 
into    the    following  sections: 

1.  PURPOSE 

This    section      describes    the      intended   purpose      of 
the    program  or   subroutine. 
2.     ass'j:iptions 

This    section    lists    any    assumptions    made    which   are 
applicable   to   the  program   cr   subroutine. 

2.  EQUATIONS 

This  section  lists  the   equations  utilizad  within 
the  main  program  cr  subroutine.    The  primary  refarances  for 
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tha      equations      nssd      are    Aerodvnaaics      and      V/STQL      ?liah':> 
Rafarence      3.1,        [jeli^opter  Thsorj,         Srf^rsncs      3.2,         a-i 
Aircraft    Performance,    Rsference    3.3. 
a.       FLOWCHART 

Both  the  handheld  computer  and  the  13^1  3033 
normally  execute  insrruct ions  in  a  program  in  a  sequential 
manner      unless   it      is      instructed      to    do      o-herwis=.  This 

section  will  graphically  represent  -he  =t=p  by  step  method 
used  to  solve  the  problem  as  well  as  th^  flow  of  con-rol 
beiween  the  various  parts  of  the  prograa.  In  a  flDwchar-, 
different  types  of  operations  are  indicated  oy  differen- 
shaoed   boxes    as    illustrated  below: 


Oval 


For    star-    or    s-oo. 


Rectangle 


For    a   oalcula-ion   or 


process    other    than    a 
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Mcdifiad 

Rectangle 


For    the    execu*:icn    cf 


a   suorcutir.  ?. 


Diamond 


For    a    decision. 


Parallslogr a.-n      For    inpa-    or   ou-pu* 


O 


Small    Circle        For    an   on    oaqe 


connrCtion    wrxe: 


flowchar-    continurs 


OQ    -ae    same    page,    or 


when  it    is  difficult 


to  connect  two  boxes 


Pentagon 


For   connection  when 
a  flowchart  continues 
to  another  page. 
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5.       PROGRAMS    AND     SUBROUTINES 


This 


1     i   =•  *■ 


5sc--iQn        lists      all        tas         prcgrans 


sabrcuT-ines    which  must,      be   prasenx    in    "cha    HP      41-C7   prior   -o 
program    axecutio?.. 

6.       PROGRAM    LISTING 

This    seo-^-ion    contains    tha    HP   4  1-CV    lisiing    of    'ha 
program    or    subroutine. 

3.        QUICK    REFERENCE    TABLES 

The    tablas   in  this   section   will   ba    usaful    to    th=    nc:^'-   '^.' 


'-i  w  tr  «       w  1. 


the      HP      U1-CV   handheld      calculator      as      a    source      of      quick 
rafaranca.  Table    1      is      an      alphabatioal   lisiing      of      all 

calculator    displays    with    an      explanation    of   thair    rrspectiva 
meaning.  Tabl?  2      lists  all      storage      ragisters   used      and 

describes    how   they  ara   utilized. 
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I-      HiiilzCV    Displays 


TABL2    I 


DISPLAY 
a    =    ? 
b    =    ? 
c    =    ? 
DA    =       ? 
Cic    =   ? 

F?A    =    ? 

h    =    ? 

aCRO=? 


R    =    ? 
TWIST    =    ? 
7?(KT)     =     ? 
VT    =    ? 
■-J    =    ? 
COLL    = 


DEFINITIOM 
Slops   of    airfoil    lift   carv9. 
Number  of    main   rotor    biaies. 
Chord   length    of    main    ro-or. 
Density    altitude. 
Wain    rotor  cosffici^n-    of    drag 
at   zero    lift. 

Equivalent   horizon-al    flat    pla- 
area    of    the   helicopter. 
Heigh-    of   the    main    rotor    above 
the    ground. 

Critical  Jiach  number  for  rhe 
coefficient  of  1  if -  equal  -o 
zero. 

Main    rotor  radius. 
Geometric    -wist    of    the    r^-or. 
Forward    velocity. 


Ret  or    tiD    soeea. 


Aircraft    aross   weiani. 


^ain    rotor   collec-ive    oitch, 


UNITS 


raa 


•I /A 


^i. 


N/A 


N/A 


rad 

Kt 

f-Vs 

lbs 
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CYCLIC  =  :iain    rotor  cyclic   pi-ch.  rai 

?I    =  Induced    power   raquirad  hp 

compensated  for    tip   losses 

and    ground  effect. 
PM    =  Power   required   due   zc  hp 

compressibility    effects. 
?0    =  Profile    power    required.  hp 

PS    =  Power   required   due    -o    -a=  hp 

retreating   blade    s-alling. 
V3E    =  Maximum    endurance   velocity.  Kt 

VF    >    v:iAX  The    forward  veloci-.y   tha-    has  N/.\ 

been    input  is    larger   thar.    the 

one    calculated   by   subroutine    VHAX. 
VMAX    =  Maximum    forward    velocity.  Xt 

VS    =  Initial    estimate    for    finding  Kt 

onset   velocity. 
XS    =  Radius   outboard    of    whicn    the  N/A 

main    rotor  is   stalled. 
X3    =  Padius   inboard  of   which   rctor  N/A 

blade   stall  may    be    presen-. 
a9  0D    =  Angle   of    attack:    of   the    advancing        deg 

blade   at        =    90    degrees. 
a270    =  Angle    of    attack    of    the   rBtr=ating      ieg 
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blade   at    i|;    =    270    degrses. 


2.      HZiilzCV    Hesistsr    Utilization 


TABLE     II 


STORAGE 

REGISTER 

00 

01 

02 

05 

06 

07 

09 

10 

1  1 

12 

13 

14 

15 

16 

17 

18 

STORED    QaANTITY/OSE 

Ground    affect  ratio    (GE)  . 

Maxiaam    forward   veloci-y    (VMAX)  . 

Stall    onset    velo::ity    (VS)  . 

Rotor    radius    (R)  . 

Number    of    blades    ( b) . 

Zero-lift    drag   coefficient    (Cio)  . 

Rotor    height  above    rha   ground    (h)  . 

Aircraft  gross   weigh-    (W)  . 

Air  density    (RHO)  . 

Lift  curve    slope    (a) . 

Rotor    lip    velocity     (VT)  . 

Coefficient    of    thrust    (CT)  . 

Tip-Loss   factor     (BTL)  . 

Main   ro-or    induced    power    (PI). 

Ro-or    height,  to    rotor    dianieter    ratio     (h/d)  . 

Compressibility    power    (PC). 
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19  Soiiii-.y    (SB). 

20  Indursd    vslccity    (VI). 

21  Main  ro-or    profile    pow^r    (P3). 

22  Advance   ratio    (MU)  . 

23  Stall    power    (PS)  . 

24  Advance   ratio  squared    (."1U2    )  . 

25  Forward   velocity    (  VF)  . 

26  Equivalent    hcriz:;n-al    flat    plate    area    (??A)  . 

27  Maximum   2-D    lift   coefficient    (CLMAX)  . 

28  Main  rotor    parasite   power    (?P)  . 

29  Main   rotor    geometric   twist     (TWIST). 

U1  Angle    of   attack    at    270    degrees     (A273)  . 

42  Angle    of   attack:    at    90    degr^is    (A90)  . 

43  Sonic    velocity    (S7EL). 

44  Critical  Mach  number    for   coefficient    of    lif- 
equal    to   zero    (MCRO)  . 

45  Density   altitude, 
OTHERS  Scratch   pad    calculations. 
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C.        FHOGHAII    DOCULISNTATION 

This    s€C-icn   contains    the    nscassary   aocumsnta- ior.    for 
the    H?  U1-CV    computer   program.      The    main    program   as   well    a: 
all    cf  the    subroutines  used  in   the    solunion   of    -ha   problem 
are    ouilinsd    in    this    section. 
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(a)  PURPOSE 

This   program    calculates    the   dynamic   parameters    of 
the    main    rotor   which    are    necessary    for   calculating   the 
total  main   rotor   power   required    in    forward,    high   speed 
straight    and   level   flight.      I-    additionally   controls   the 
execution    sequence   of   the    various   subroutines    which   are 
used   to    calculate  main   rotor   power    required,    in   tens 
of    horsepower,    as   well    as    maximum    forward    velocity,    stall 
onset  velocity  and  velocity  for    best   endurance   in   kno-s. 

(b)  ASSUMPTIONS 

(1)  All    angles    are    small. 

(2)  Steady    flow   through   the   rotor. 

(3)  All   rotor    blades  are   rectangular    (non-tapered) 
with  only    uniform  twist    being   possible. 

(U)       Only   the   first    hamonic    of   flapping    is   necessary 

for    calculating   power    required. 
(U)       The    effective    dimensionless    radius    can    be 

apprcxima t edy   by   the   tip-loss  factor. 

(6)  The    thrust    vector  passes    through    the    C.G,. 

(7)  The    static    stall  angle    for   blades    on    most 

helicopters   can    be    approximated    by    the    angle    at 

which    CLl^AX    occurs    for    the   2-D    airfo-l. 
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(8)       Th=      6    3   affsct,   or    the   result    of    cccicir.g   th' 
flapping   axis    of  the   blada   sc   t  ha-    its    pitch 
varies    as    the    blade    flaps   is    zaro    (i.e.,    no 
lateral   flapping  is    present). 


(C)        EQUATIONS 


u    =       Vp/V^ 

V     (ft/s)    =    V    *    1.68394 

F 
T1    =    .5(B^2     ^    .5    .12) 

T2    =    .3348    3    +    .  5   ^23 

TL  TL 

T3    =    .25B„2C  32     +  ^^2) 
TL      TL 

TU    =    .5y  (B     2    +    ,  25   u^^ 

TL 
D1    =     (B    2    -    .5    y2) 

TL 
All    =    4  (.5    B    2    -    -3    /8)  /D1 
TL 

A12    =    BuB      /3D1 

A13    =    2yB   2     /    D1 

A14    =     (B    2     +     1.5    y2)     /     D1 

TL 

(2CT/aa)    =     ATI    *    9   r 2    +     9„ T3    +    9    T4 
^  '  0  T  2 

0    =    XA11    <•    9    A 12    +9„A13    +    9    A14 

0  T  2 

Amax    =    CLMAX/CLA 

a270    =    9         -9       +9+(A/1+y) 
0  2  T 

r       =         AMAX    -0      -0 
0  2 

Cg   =        yr   +    X 
33  =    -  u9^-   r 

^    =    -^S-       ^^S"     -   ^^T-S)*"    /    29^ 
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d.   FLOWCHART 


Start 


fnput : 


^n. 


subroutine 


DEN 


suDrout Tne 


DW 


subroutine 


k 

PIT 

subroutine 

1 

DOl 

subrout ine 
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2R1_ 


sudroutTne 


jvi. 


sudroutfne 


subroutine 


set  flag 
02 


DISPLAY: 
VF  >  VMAX 


, eg  ■ 

sudroutTne 


GO  TO  A 


CT 


subroutine 


calculate: 
T1,T2,T3, 

T4 
1 


calculate; 

A11,A12, 

A13,Ai4 


LAMB 


subroutine 


Jl 


calculate; 

xs.xo 


CPS    =    0 
(VF    <    VS) 


a90 


subrouti ne 


XRL 


subrouT  f  ne 
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NO 

=    0 

1 

r 

CG 

subroutfne 

r 

GO    TO   A 

e.   PROGRAMS  AND  SUBROUTINES  USED 

♦♦  CT  ♦♦ 


n    pjj  n 
n    poi  ♦• 

tf  V5  n 
"  VMAX  ' 


"  LAMB  ' 
♦•  CPS  - 
♦•  CPC  " 
*•  a90  ♦• 
-  CNG  *• 
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f.   PROGRAM  LISTING 


01+LBL  "i4BS- 
92  'R—'" 

83  PROHPT 

84  STO  85 

85  •b=''" 

86  PROHPT 
97  STO  86 

88  •C='^- 

89  PROHPT 
19  STO  84 

11  •r:dO=''' 

12  PROHPT 

13  STO  87 

14  ■H=^'- 

15  PROHPT 

16  STO  18 

if     ?  !  —  ■ 

18  PROHPT 

19  STO  13 
20  ■VF<KT>=?" 

21  PROHPT 

22  1,68894 

23  * 

24  STO  25 

25  •pppr'?- 

26  PROHPT 

27  STO  26 

28  'H='" 

29  PROHPT 
38  STO  89 

31  "5=^- 

32  PROHPT 

33  STO  12 
34  •CLHfl>;='?- 

35  PROHPT 

36  STO  27 
37  ■THIST=?- 

38  PROHPT 

39  STO  29 
48  •HCR8='- 

41  PROHPT 

42  STO  44 


43  ■Dfi=^- 

44  PROHPT 

45  STO  45 


4€*LBL  'RGH" 

47  RCL  25 

48  RCL  13 

49  / 

58  STO  22 

51  Xt2 

52  STO  24 

53  XtQ  •SD- 
54  XEQ  -DEH- 

55  FS'5  83 


57  KEQ 

58  XEQ 


■PIT' 
■POl" 


59*L8L  95 
68  ;0  'PPl' 

61  flIiV 

62  FS"?  93 

63  GTG  86 
fA  'm   -VBE" 

65  XEQ  -VHflX- 
66  mDV 

67*LBL  86 

68  KEQ  "CT- 

69  RCL  15 
79  Xt2 

71  STO  39 

72  EHTERt 

73  RCL  24 

74  2 

■■J  ' 

76  ^ 

77  STO  31 

78  2 

79  / 

88  STO  33 
81  RCL  15 


82  3 

33  nx 

84  3 

85  / 

86  EHTERt 

87  RCL  24 
38  RCL  15 
89  * 

98  2 

91  / 

92  •*■ 

93  STO  34 

94  RCL  39 

95  RCL  24 

96  + 

97  RCL  39 

98  * 

99  4 
199  ■' 

191  STO  35 

182  RCL  39 

183  RCL  24 

194  4 

195  / 
186  + 

107  RCL  22 
188  * 
199  2 
119  / 

111  STO  36 

112  RCL  22 

113  RCL  30 

114  * 

115  2 
lib  / 

117  RCL  22 

<  1  0  1 

119  YtX 
128  8 
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Ibb  RCL  12  287  - 

124  *  ^^^  ■''  288  / 

125  RCL  38  168  RCL  19  299  STO  32 

126  EHTEFt  169  /  218  F!?'?  87 

127  RCL  24        ■  178  STO  42  211  QlQ   «7 
123  2  171  XEQ  -LfiHB-  212  R-S 
129  /  i72  RCL  83  21 3^" CYCLIC^ 
138  -  173  RCL  33  214  fiRCL  X 

131  STO  32  174  *  215  WltU 

132  RCL  38  175  CHS  216  Sinp 

133  *  176  RCL  42 

134  /  177  *■  217*LBL  87 

135  STO  37  173  RCL  29  218  RCL  31 

136  RCL  22  179  RCL  35  219  EHTERt 

137  RCL  15  188  *  228  RCL  36 

138  *  131  -  221  RCL  32 

139  8  182  STO  31  222  * 
148  *  183  RCL  83  223  - 

141  3  184  CHS  224  RCL  34 

142  /•  185  RCL  37  225  -^ 

143  RCL  32  186  *  226  STO  33 

144  /  187  RCL  29  227  FS?  83 

145  STO  38  183  RCL  39  228  GTO  88 

146  RCL  22  189  «  229  R-B 

147  2  198  -  238  •caLL=" 

148  *  191  STO  32  231  PRCL  X 

149  RCL  38  192  RCL  34  232  PVIEH 
158  *  193  *  233  STOP 

151  RCL  32  194  ENTERt  234  PDV 

152  /  195  RCL  31 

153  STO  39  196  RCL  38  235*LBL  88 

154  RCL  38  197  *  236  RCL  33 

155  EHTERt  198  -  237  RCL  32 

156  RCL  24  199  RCL  34  238  - 

157  1.5  288  EHTERt  239  RCL  29 

158  *  281  RCL  48  248  * 

159  +  282  *  241  EHTERt 
168  RCL  32  283  RCL  36  242  RCL  83 

161  /  284  EHTERt  243  1 

162  STO  48  285  RCL  38  244  RCL 

163  RCL  14  286  *  245  + 

164  2 

165  * 


'I'i 
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246 

/ 

247 

+ 

248  STO  41 

249 

FS?  1 

33 

25e 

GTO  ! 

89 

251 

R-n 

252 

•<i27' 

R=- 

253 

PRCL 

X 

254  fiVIEy 

255 

STOP 

256 

anv 

257  KtQ  • 

m 

258  flDV 

259 

RCL 

27 

266  RCL 

12 

261 

,-•• 

262 

RCL 

33 

263 

- 

264  RCL 

32 

265 

+ 

266 

STO 

3fi 

267 

RCL 

22 

268  CHS 

269  RCL 

29 

27fi 

« 

271 

RCL 

30 

272 

- 

273 

STO 

31 

274 

RCL 

38 

275 

RCL 

276 

« 

277 

RCL 

83 

278 

•f 

279 

STO 

35 

288 

4 

281 

CHS 

282 

* 

283 

RCL 

29 

284 

* 

285 

STO 

•J -J 

286 

RCL 

31 

287  Xt2 

288  STO  36 

289  ENTERt 
298  RCL  35 

291  CHS 

292  JOV^ 

293  SF  81 

294  FS^  81 

295  GTG  83 

296  RCL  36 

297  RCL  35 

298  + 

299  SSRT 
380  STO  38 
301  RCL  31 

382  CHS 

383  + 

304  RCL  29 

305  / 

306  2 

307  / 

388  STO  38 
309  ■XS=' 
318  RRCL  ?; 

311  flVIEy 

312  STOP 

313  RCL  31 

314  CHS 

315  RCL  38 

316  - 

317  RCL  29 

318  / 

319  2 

320  / 

321  STO  88 

322  "X0=' 

323  flRCL  X 

324  flVIEH 

325  STOP 

326  PIIV 


327*LBL  83 
328  XEQ  'CPS' 

329  RDV 

330  CF  81 

331  CF  82 

332  CF  83 
333  XEQ  -CPC- 

334  9BV 

335  Prv 

336  XEQ  -PTl- 

337  KE9  -CHG- 

338*LBL  89 

339  XEG  -VS- 

340  END 
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2 .      Solidj-t-Z 


a.  PURPOSE 

Thfs  subrouttne  calculates  the  ratto  of  total  blade  area 
to  the  total  rotor  atsc  area. 

b,  ASSUMPTIONS 
None 


c.  EQUATIONS. 


b  ♦  c  *  R 


b  *  c 


PI  ♦  R 


PI  *  R 


d.   FLOWCHART 


start 


calculate 
solfdtty 


6U 


»i.      PftOOJAMS    AND    5UBP0UTINES    USED 

*    WBS    ♦ 
f.      PROGftAM   LISTING 


31*LBL  "SB- 
Be  RCL  Mb 

83  RCL  9* 

84  * 

85  RCL  95 

86  ■' 

87  PI 
08  ■' 

99  STO  19 
19  EHT. 


bD 


3.      Down wash    Velocity 

a.     PURPOSE 

The  purpose  of  thfs  subroutine  fs  to  compute  the  Induced 
velocfty  of  a  rotor  system. 

b.  ASSUMPTIONS 

Steady  flow  through  the  rotor  system. 

c.  EQUATIONS 

U 

w 


2PA,  V 


d.   FLOWCHART 


start 


calculate 
downwash 
velocity 


store  DW 
tn  R32 


return 

e.   PROGRAMS  AND  SUBROUTINES  USED 
•  WBS** 
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f.      PROGRAM  LISTING 


91* 

LBL 

„y. 

82 

RCL  ; 

25 

fi7 

X=8^ 

84 

GTO 

85 

85 

RCL 

18 

9b 

0 

8? 

RCL 

11 

88 

* 

89 

RCL 

25 

18 

* 

11 

PI 

12 

* 

13 

RCL 

95 

14 

At2 

15 

* 

16 

/ 

1? 

STQ 

64 

18 

GTO 

8b 

15 

♦LBL 

85 

28 

1  SF 

93 

21 

.♦LBL 

86 

22  EHI 
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^«      Coefficient    o_f  Thrust 

a.     PURPOSE 

Thfs  subroutTne  calculates  the  coefftcfent  of  thrust  for  a 
arbitrary  rotor. 

d.   ASSUMPTIONS 

Steady  flow  through  the  rotor  system. 

c.   EQUATIONS 


CI     =   '*' 2 

^'    pA  V  ' 

^  0  T 


d.   FLOWCHART 


Start 


calculate 

coeff  of 

thrust 


Store  Tn 
Ria 


return 
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e.      PROGRAMS    AND    SUBPOUTINES   USED 

•    WB5    ♦ 
f  .      PHOGftAM   LISTING 


ei*LBL  -C" 
82  RCL  95 
93  XtZ 
64  PI 
95  * 

86  RCL  11 

87  * 

88  RCL  13 

89  Xt2 
18  * 

1  \     1  .-v 
i  1  i  -■  n 

12  RCL  19 

13  * 

14  STO  14 

15  EHB 
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5«      Induced    Power 

(a)  PURPOSE 

This    subroatine    calculates    the    pDwer    r9quir==i    by   th'?    ~oi' 
-0    produce   thrust  at. hover  and   forward    flight.       Addi- ionally , 
-his    subroutine    corrects    for   tip    losses    (losses    in    lift    at    the 
tips    due    tip   votices)    as    well   as    ground   effect. 

(b)  ASSUMPTIONS 

Steady    flow   through    the   rotor    systeai. 
(C)        EQUATIONS 


PI 


2v| 


'V£Y       -     1 


2vJ 


1.0     -V2    CT 


GE      = 


PIT   = 


-0.  1276(h/D)  ♦+O708  (h/D)  3 
-1.4  5  69(h/D)  2+1  34  3  2  (h/D)  +0.514  7 
(1/B       )*(GE)  *PI 
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d.   FLOWCHART 


c 


start 


calculate 

tip  loss 

factor 


YES 


NO 


G.E.  =  1.0 


calculate 
GE   ratio 


calculate 

Induced 

power 


compensate 

for  tip 

losses 


compensate 
fo  G.E. 


return 
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e.  PftOO<AMS    AND    SUBPOUTINES   USED 

•    WB5    • 

f.  PHOGftAM   LISTING 

33  .788 


81n 

.BL  ■ 

82 

RCL 

14 

83 

y 

84 

* 

85 

SQR^ 

86 

RCL 

86 

87 

/ 

88 

CHS 

89 

1 

18 

* 

11 

STO 

15 

12 

RCL 

89 

RCL 

85 

14 

-J 

15 

* 

16 

/ 

17 

i.  - 

STO 

4  •? 

18 

1.5 

19 

XOY 

28 

x>y^ 

21 

GTG 

85 

22 

RCL 

17 

23 

1.3432 

24 

X 

'?5 

RCL 

1  J 

26 

"t2 

27 

-1.4569 

28 

* 

29 

■•• 

38 

RCL 

17 

31 

■1- 

3? 

T+V 

34 

* 

35 

■t 

36 

RCL 

17 

37 

4 

38 

'in 

39 

-.1276 

48 

* 

41 

t 

42 

.51' 

^7 

43 

4- 

44 

STO 

88 

45 

GTO 

86 

46<eL 

85 

47 

i 

48 

STO 

88 

49*LBL 

86 

58 

RCL 

25 

51 

?;t2 

52 

STO 

32 

53 

RCL 

18 

54 

^, 

55 

/ 

56 

RCL 

11 

57 

/ 

58 

PI 

/ 

68 

RCL 

85 

61 

st2 

62 

/ 

63  SORT 


64  STO 

28 

65  Xt2 

66  1/X 

67  RCL 

32 

68  * 

69  2 

78  / 

71  CHS 

72  STO 

34 

73  ^t? 

74  1 


76 

SQRT 

77 

RCL 

34 

78 

■<■ 

79 

SQRT 

88 

RCL 

18 

81 

* 

82 

RCL 

28 

83 

* 

84 

558 

85 

86 

RCL 

87 

* 

88 

RCL 

88 

89 

* 

98 

STO 

16 

Q1 

"PI= 

■ 

92 

ORCL 

93 

qviEH 

94 

STOP 

95  END 
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6.  2l.2^L2i^    poy  ar 

a.  PURPOSE 

This  subrouttne  calculates  the  oroffle  power  required  for 
forward,  straight  and  level  flight  In  terms  of  horsepower. 

b.  ASSUMPTIONS 

Steady  flow  through  the  rotor  sustem. 


C.   EQUATIONS. 


P   -   J  Cdo  p  A  vMl  ♦  4.25  u^) 

0  0  T 


d.   FLOWCHART 


4400 


(hp) 


start 


^ 


calculate 

profile 

po\yer 


store  PO 
In  R21 


return 


e,   PROGRAMS  AND  SUBROUTINE  USED 
*  W8S  " 


73 


f.      PROGRAM  LISTING 


01*L 

.BL  -FOl 

92 

RCL  19 

83 

RCL  87 

8* 

* 

85 

RCL  11 

% 

* 

87 

RCL  85 

88 

1*7 

89 

* 

18 

PI 

i  i 

* 

12 

RCL  13 

13 

3 

14 

'in 

15 

* 

16 

8 

77 

/ 

18 

•jiU  Hi. 

1Q 

A  -" 

RCL  24 

28 

4.25 

?' 

wl 

* 

22 

1 

23 

+ 

24 

RCL  21 

25 

* 

26 

558 

27 

/ 

28 

STG  21 

29 

■pri=- 

38  flRCL  K 

31 

aviEy 

32 

STOP 

33 

tm 

7U 


7-      ZH^si_t9    Power 

a.  PUPOSE 

ThTs  subroutfne  calculates  the  parasite  power  requfred  tn 
forwar(3»  straight  and  level  flight. 

b.  ASSUMPTIONS 

Steady  flow  through  the  rotor  system. 

c.  EQUATIONS. 


PP»  Jll 


1100 


(hp) 


d.   FLOWCHART 


/"  start 


calculate 

parasite 

power 


store  PP 
In  R28 


return 


§.  PROGRAMS  AND  SUBROUTINES  USED 


WBS 
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f.      PROGRAM   LISTING 


81*1 

BL  "P^l 

82 

RCL  11 

93 

RCL  26 

84 

* 

85 

p 

86 

* 

8? 

RCL  25 

88 

3 

89 

TtX 

18 

* 

11 

558 

12 

/ 

13 

STO  28 

14 

■PF=* 

15 

flRCL  X 

16 

hviem 

17 

STOP 

18 

END 

76 


Q«   ^.axinum  Fee  ward  7  eloci-]^ 


a.  Purpose 

This  subrouttne  calculates  the  power-l Tmfted  maxtmum  speed 
of  the  spec tf ted  he It copter, 

b.  ASSUMPTIONS 

a.  The  power- 1  tint ted  maxtmum  veloctty  may  be  esttmated 
by  neglecttng  the  vartattn  of  tnduced  power  and  prof  tie 
power  wtth  speed. 

b.  Power  required  to  hover  ts  approxtmately  equal  to 
power  requtred  fo  maxtmum  speed. 

c.  Steady  flow  through  the  rotor  system. 


c.  EQUATIONS. 


VMAX  =    V,  [_!_  ]  ^^^ 


W 

V         -      r     3  1/2 

*        '■  2  p  A  -^ 

0 


77 


D.   FLOWCHART 


start 


calculate 

tnduced 

velocHy 


calculate 
max  fwd 
velocHy 


NO 


return 


YES 


set  flag 
02 


NO 


YES 


return 


dfsplay 
VF  >  VMAX 


execute 

subrouttne 

change 


78 


PROGRAMS   AND    SUflPOU TINES   USED 

•    W85    ♦ 

PHOGftAM   LISTING 

9i*LBL  "VHfiX' 
82  RCL  2b 
93  PI 

85  RCL  85 

86  ?:t2 

87  -^ 

88  l/K 

89  4 
18  * 

11  .3333334 

12  YtK 

13  RCL  28 
1*  * 

15  STO  91 
lb  1.68S4 

IS  •VHfiX=- 
1?  QRCL  X 
28  flVIEy 

21  STOP 

22  RCL  81 

23  RCL  25 

24  Km 

25  GTO  12 
2b  X<=Y^ 
27  GTO  13 

28*LBL  12 

38  aVIEM 
31  aEQ  ■CNG- 

32*LBL  13 
33  EHD 


79 


9-      Best    Er.  dun  nee   V  excci  tv 

A,  PURPOSE 

To  calculate  the  value  of  velocfty  correspondTng  to 
mfnfmum  power  (t.e.,  best  endurance  velocTty  and/or  best 
rate  of  c1 fmb) , 

B,  ASSUMPTIONS 

1.  Steady  flow  through  the  rotor  system. 

2.  The  varfatton  of  prof  fie  power  wfth  forward  velocfty 

fs  negltgfble. 

C,  EQUATIONS 


,.  r      U .      / A s1/2t    1/2 

"^'^     '    ^    p  *  A      *  ^3   *   f    ^      ^ 


(ft/s) 


D.      FLOWCHART 


start 


calculate 
VBE 


convert 
VBE  to 
knots 


return 


30 


e.  PftOGHAMS    AND    SUBPOUTINES   USED 

♦    WBS    ♦ 

f.  PROGRAM   LISTING 


81*1 

SL  =VBE" 

Ar 

RCL  85 

■:^7 

xt2 

iS4 

PI 

85 

» 

96 

3 

87 

/ 

88 

RCL  26 

89 

/ 

18 

SQRT 

11 

RCL  18 

1-? 

* 

RCL  11 

14 

/ 

15 

RCL  85 

16 

'm 

1-? 

/ 

18 

PI 

19 

/ 

28 

SQRT 

21 

1.68894 

22 

/ 

23 

•VBEK=- 

24 

qPCL  X 

25 

RVIEil 

26 

STOP 

27 

EHB 
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10.       S;tall    Onggt    Velocity 


a .  PURPOSE 

Thts  subroutTne  gfves  the  user  of  the  HP  41-CV  an  tnftfal 
approxtmatfon  for  the  velocUy  at  which  the  retreatTng  blade 
angle  of  attack  fs  approximately  equal  to  the  Is  approximately 
equal  to  the  static  stall  angle  of  the  rotor  blade. 

b.  ASSUMPTIONS 

1.  Steady  flow  through  the  rotor  system. 

2.  Stall  onset  velocity  Is  approximately  equal  to  the 
velocity  for  best  range  (I.e.,  minimum  P/V), 

c.  EQUATIONS 

p   =   (-4.25  a Cdo  A^V^  )/(4f) 

r   -   (-a  Cdo  A  vj)/(4f) 


a 


=  -1/3  *  [^ 


b  =   1/27  *  (2  ♦  p^+  27  *  r) 


A  »  C  -b/2  ^  (bV2  *   aV27)^^^]  ^^^ 
B  »  C  -b/2  -  (bV2  *  i  /21)^'^1  '^^ 
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d.   FLOWCHART 


start 


calculate 
A  and  B 


fTnd 

largest 

cuDfc  root 

1 

' 

VSl  =•  A  ♦ 
8 

' 

VF  -  VSl 

<D 


calculate 
a270 


YES 


VF  =  VF  ♦ 
5 


VF  =  VF-5 


e.   PROGRAMS  AND  SUBROUTINES  USED 


WBS 


YES 


VS  -  VF 


return 


83 


f.      PROGRAM   LISTING 


81*LBL  "VS" 
62  FS''  93 
03  GTQ  10 

84  RCL  19 

85  RCL  87 

86  * 

88  * 

89  RCL  85 

18  xn 

11  * 


12  4 

4  7 

I  J    ■■ 

li  RCL 

26 

15  / 

lb  STO 

83 

17  RCL 

<7 

18  * 

19  4.25 

28  * 

21  CHS 

22  STO 

81 

23  RCL 

i  7 

i-' 

24  3 

25  YtX 

26  RCL 

83 

27  * 

28  CHS 

29  STO 

82 

39  27 

32  RCL 

81 

33  3 

34  VtK 

■h  * 


38 

27 

39 

/ 

48 

STO 

31 

41 

RCL 

81 

42 

Kt2 

43 

^A 

1/X 

45 

* 

46 

STG 

38 

47 

RCL 

31 

48 

Kt2 

49 

4 

58 

/ 

51 

RCL 

30 

52 

7 

53 

■ftx 

54 

27 

55 

/ 

56 

■i- 

57 

SORT 

58 

STO 

33 

59 

RCL 

31 

68 

CHS 

61 

.- 

62 

/ 

63 

+ 

64 

J3J 

334 

65 

YtX 

66 

STO 

32 

67 

RCL 

31 

68 

CHS 

69 

2 

78 

/ 

71 

RCL 

33 

72 

- 

^•7 
5'  -J 

.333334 

74 

TtK 

75 

STO 

34 

76 

RCL 

32 

77 

4- 

78 

STO 

82 

79 

.25 

88 

* 

81 

RCL 

82 

82 

+ 

83 

STO 

25 

84 

SF  1 

53 

85 

GTO  ' 

■qCH 

86*LBL 

18 

37 

RCL 

'■■-? 

if 

88  RCL 

12 

89 

/ 

98 

STO 

31 

91 

RCL 

41 

9? 

_ 

Q7 

'i^? 

94 

SQRT 

95 

.881 

96 

'i<yf 

97 

K<=Y'? 

98 

GTO  11 

99 

RCL  31 

188  RCL  41 

181 

m'> 

102 

GTO  12 

183 

RCL  41 

184 

RCL  31 

185 

X)^? 

186 

GTO  13 

107*LBL  12 

188 

RCL  25 

189 

c; 

118 

- 

ill 

STO  25 

112  GTO  -PGH' 

113*LBL  13 

114 

RCL  25 

115 

5 

116 

-i- 

117 

STO  25 

118  ! 

^TO  -OGN' 

119*LBL  11 

126 

CF  93 

121 

RCL  41 

122 

R-Ii 

123  ' 

■iSTOLL= 

124 

PRCL  X 

125 

fiVIEH 

126 

STOP 

4  '\-; 

RCL  25 

128 

1.68894 

129 

/ 

138 

■VS=" 

131 

QRCL  X 

132 

PVIEH 

1  "7"' 

STOP 

134 

ENI! 

84 


11.      Injlo^i   P-atio 

(a)  PURPOSE 

This    subroutins    calculates   the   ratio    of   -cha   r.3t    veloci-.y    u? 
through    the   rotor  systsni    to  tha    tip   spaai. 

(b)  ASSUMPTIONS 

Sready    flew   through    the   rotor   sys-r.-n. 


(C)        EQUATIONS 


X      =  -    /CT/2 

X      -  '    ,  CT  +    y   Tan  a 

Dp   =  (PP*550)/Vp 

a  2   =  -Tan     {C  /W) 


85 


d,   FLOWCHART 


/ 


/  Start 


y 


yes 


return 


calculate 
1  m  1 1.  a  1 


no 


YES 


Y 

return 


calculate 
a3 


calculate 

X 


NO 
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e.  PftOGftAMS    AND    SUflPOUTlNES   USED 

•    WBS    • 

f.  PROGftAM   LISTING 

81 +181  "LfiHB" 


84  / 

85  S8RT 
96  CHS 

87  STO  93 

88  =1 

89  RCL  22 
lg  x<=V? 

11  GTO  86 

12  RCL  28 

13  558 

14  * 

15  RCL  25 

16  /■ 

17  RCL  18 

19  STO  23 


39  RCL  83 
48  - 

41  Xt2 

42  SQRT 

43  .88885 

44  K>Y? 

45  GTO  82 

46  RCL  88 

47  STO  83 

48  QTG  81 

49*LeL  82 
58  RCL  88 

51  STO  83 

52*LBL  86 
53  EHD 


2e*LBL 

81 

21 

RCL 

14 

-t'^ 
c^ 

CHS 

?7 

RCL 

83 

24 

Xt2 

25 

RCL 

24 

2  b 

■i- 

27 

SQRT 

28 

2 

29 

* 

38 

/ 

•71 

EHTERt 

32 

RCL 

"■7 

33 

R-D 

34 

TPH 

-71: 

•_»■-■ 

RCL 

00 

36 

* 

37 

- 

38 

STO 

88 

87 


12.       Anal*   2l  Kllk^k   ai  11  DeQ-^^3 


a.  PURPOSE 

Ihfs  subroutine  calculates  the  angle  of  attack  at  the 
azfmuthal  posftTon  of  90  degrees. 

b.  ASSUMPTIONS 

a.  Steady  flow  through  the  rotor  system. 

b.  Blade  oscTllatTons  are  perTodtc  tn  nature. 

c.  Only  first  harmonfcs  of  flapping  are  necessary  for 
calculating  angle  of  attack. 

d.  The  thrust  vector  passes  through  the  C.G. 

e.  Only  uniform  twist  of  the  rotor  blade  Is  possible. 

c.  EQUATIONS 

X 


a90   =  e  +  e  +  e  + 


0    2     T   (  1  +  y  ) 


88 


d.   FLOWCHART 


start 


recal 1 

col lectTve 

anqle 


recal 1 
cycl Tc 
angle 


recal 1 
1 fnear 
twfst 


calculate 
a90 


Store  Tn 
R42 


return 


39 


e.  PROGRAMS  AND  SUBROUTINES  USED 

♦  WBS  * 
"  a90  •• 

f.  PROGRAM  LISTING 


Sin 

.BL  -m- 

82 

RCL  33 

93 

RCL  32 

84 

+ 

85 

RCL  29 

86 

+ 

87 

ENTERt 

88 

RCL  83 

89 

18 

RCL  22 

i  i 

+ 

12 

/ 

13 

+ 

14 

STQ  42 

R-D 

16 

"i98I>=° 

i  "7 

1  l' 

wRCL  X 

15i 

fiVIEM 

19 

STOP 

28  END 

90 


13.   Cgmprggsibilit y  Powjr 

a .  PURPOSE 

This  subrouttne  calculates  the  power  requfred  due  to  com_ 
pressfbtlfty  on  the  mafn  rotor  sustem  tn  forward .straTght 
and  level  flTght  tn  terms  of  horsepower, 

2.   ASSUMPTIONS 

a.  Steady  flow  through  the  rotor  system. 

D.  The  compresstbTl Tty  losses  can  be  expressed  as  a  function 
of  the  amount  by  whfch  the  drag  divergence  Mach  nunber  Ts 
exceeded  at  the  tip  of  the  advancing  blade. 

c.   EQUATIONS 


^to      a 


1 

o 


AM    -    M   -  M  ,.  0.06 
d         t      crft 


M  .   -   0.71  -  2.3  ♦  a90 
erft 


C     -     *  [0.012  *  AM   >  0.1  *  (AM^)  *] 
PC  d  d 


P   •   C   *p*PI*R'*V*      (HP) 

H         PC    "^  t 

550 


91 


d,   FLOWCHART 


start 


calculate 

t!p  Mach 

number 


calculate 
AMD 

> 

' 

recall  a90 

1 

' 

calculate 

crTtfcal 

Mach  No 

' 

1 

calculate 

compressfbTl tty 

power 

92 


f.      PROGRAM   LISTING 


ei*LBL  ' 

"CFC 

92 

RCL 

17 

83 

RCL 

■JP 

S4 

* 

85 

RCL 

43 

96 

/ 

97 

STO 

71 

98 

RCL 

42 

99 

RCL 

* 

n 

.u: 

.' 

12 

* 

17 

CHS 

I  \j 

14 

RCL 

44 

15 

* 

16 

STO 

32 

17 

CHS 

18 

RCL 

31 

19 

•f 

29 

.96 

"51 

- 

22 

STO 

33 

"ii 
i.-- 

9 

24 

XOt 

25 

X<='? 

'0 

26 

GTO 

99 

27 

RCL 

33 

00 

7 

v-y 

"- 

29 

TtX 

39 

i 

X  1 

31 

* 

■jC 

EHTERt 

■jO 

RCL 

33 

34 

.912 

75 

* 

37 

RCL 

19 

38 

* 

39 

STO 

34 

49 

RCL 

U 

41 

* 

42 

RCL 

95 

43 

Kt2 

44 

* 

45 

PT 

46 

* 

47 

EHTt 

■Rt 

48 

RCL 

13 

49 

/. 

59 

i  '.  - 

51 

* 

CI 

559 

/ 

54 

STO 

18 

55 

GTO 

19 

56*LBL 

99 

57 

9 

58 

STO 

18 

59*LBL 

19 

69 

■PC= 

a 

61 

flRCL 

t 

62 

PVItiii 

63 

STOP 

64 

ENI' 

93 


^^'       stall   Power 

a.  PURPOSE 

Thfs  suDrouttne  estTmates  the  addlttonal  powrer  requtred  tn 
forward,  straight  and  level  fitght  due  to  retreating  blade 
stall.   AddH Tonally,  thfs  subroutine  calculates  a  stall 
correction  factor,  k  that  corrects  for  the  special  case. 

s 

of  inboard  stal 1 tng. 

b.  ASSUMPTIONS 

a.  Steady  flow  through  the  rotor  system. 

b.  The  section  drag  coefficient  at  stall  Jumps  approxi- 
mately 0.08  at  stall  onset. 

c.  The  stalled  area  is  symmetric  about  the  270  degree 
azimuthal  position. 

d.  For  all  airfoils  considered,  the  static  stall  angle  is 
approximately  12.5  degrees. 


EQUATIONS 

-B 


2  *  G 

T 


1.0 


where. 


B  «  -u  *  o  -  r 

r  -  0.218166  -  e  +  e_ 


k  -.-( 


B  /2  *  e  +  X 

•  T     • 


1  -  X 


*  (1  -  u  )^*  (1  -  X)  *  (1  -  X^  ) 


C   =  L  24  *  pi   '"    ^    '         '"     .         • 


ff 


p  =c  ♦p*pi*r'*v' 

s     o»  T 


94 


d.   FLOWCHART 


start 


t 

NO 

<^^2eV> 

YES 

r 

w 

k  »  1.0 

s 

calculate 

stall 
correct  Ton 

' 

• 

calculate 
stall 
power 

\ 

f 

apply 

stall 

correct Ton 

<! 

1 

return 


e.   PROGRAMS  AND  SUBROUTINES  USED 

"  WBS  '• 

•  SO  n 

••  CT  •• 

»  y  It 

*  LAMB  •• 
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f.   PROGRAM  LISTING 


Ml ♦LBL  ■ 

■  rp 

02 

FS^ 

81 

83 

GTO 

96 

84 

RCL 

38 

8j 

RCL 

88 

86 

•f 

87 

^ 

88 

/ 

mQ 

1 

1 

18 

X<>' 

r' 

11 

S<=V^ 

12 

GTO 

81 

17 

i 

1 

14 

STO 

23 

15 

GTO 

82 

16*LBL 

81 

17 

i.  1 

RCL 

31 

29  RCL  29 
11 

22  RCL  38 

24  CHS 

25  EHTERt 


31 

♦LBL  82 

32 

"KS=- 

33 

QRCL  K 

34 

mtu 

7C 

•_•■-■ 

STOP 

36 

SF  82 

yy 

RCL  19 

38 

24 

39 

/ 

40 

PI 

41 

/ 

42 

EHTERt 

43 

i 
i 

44 

45 

- 

46 

Xt2 

4? 

* 

48 

1 

49 

RCL  38 

58 

- 

* 

e;^ 

i 

53 

RCL  38 

54 

At2 

'^^ 

_ 

(d 

* 

73 

550 

74 

/ 

7c: 

STO 

76  FS?  82 

77  GTO  97 

7S»LBL  86 

79  8 

88  STO  23 


81*LBL  97 
32  "PS=° 

83  flRCL  X 

84  aviEy 

85  STOP 

86  END 


56  SORT 


27  RCL  38 

■10  _ 


38  STO  23 


58  RCL 

59  * 

68  STO 

23 

61  RCL 

13 

62  3 

63  YtK 

64  RCL 

-.7 
a -J 

65  * 

66  PI 

67  * 

68  RCL 

85 

69  Xt2 

78  * 

71  RCL 

11 
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15.       Total    Power 

,      PURPOSE 

This  subrouttne  calculates  the  total  Dower  required  !n 
forward,  stratght  and  level  flfght,  to  include  stall  and 
compress t bf 1 fty  power  fn  terms  of  hosepower, 

.   ASSUMPTIONS 

Power  losses,  such  as  transmfssfon  and  cooling,  can  be 
Ignored. 


c.  EQUATIONS. 


d.   FLOWCHART 


P   =P   +P   +P  +P  +p 

T      f       0       p      9       M 


Start 


reca  1 1  P, 
requTrea 


recall  P 
requfrerf 


recall  P 
requTrecf 


recal  1. 
stall 
power 


© 


recall  P 
power  " 


sum  a  1 1 

power 

required 


return 


97 


e.  PROi;^<AMS  AND  SUBPOU FINES  USED 

•  WBS  » 

f.  PROGRAM  LISTING 


91*LBL  "PT 

62 

RCL  28 

83 

RCL  It 

84 

+ 

85 

RCL  21 

86 

+ 

87 

RCL  23 

88 

+ 

99 

RCL  18 

18 

+ 

11 

•PT=" 

12 

fiRCL  K 

13 

PVIE« 

14 

STOP 

15 

END 
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16.      pgnHJ'y/Sonic   V9 IccitY 


PURPOSE 

To  generate  values  of  atr  densTty  and  sonTc  velocity. 

ASSUMPTIONS 

(1)  Geopotentfal  altttude  (H),  and  geometrtc  altitudes  are 
equal  below  20.000  feet  (actual  ah  =  29  ft). 

(2)  In  the  troposphere  the  standard  temperature  lapse 
rate  ts  -3.57  "F  per  1000  feet. 

EQUATIONS 

e  =   T/T    »  (1.0  -  6.8753  E-06  *  H  ) 

SSL 

SVEL  =  SVEL  *   SORT(e) 

RHO  =  0.0023769  ♦  (1.0  +  HTH  •  (-.02875  *   0.000275  *  HTH) 
where, 

HTH  =  DA/1000 
FLOWCHART 


start 


calculate 

sonfc 
velocity 


calculate 

atr 

density 


return 
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e.  PftOGHAMS    AND    5UBP0UTINES    USED 

f.  PROGftAM   LISTING 


0Hl 

BL  -DEN- 

92 

RCL  45 

93  A. 

875  E-96 

94 

* 

95 

CHS 

96 

i 

97 

+ 

98 

SQRT 

99 

1116.89 

19 

* 

-1  4 

I  i 

STO  43 

12 

RCL  45 

13 

1999 

H 

/ 

15 

STO  33 

16 

.999275 

17 

* 

18 

-.92875 

19 

+ 

28 

RCL  33 

21 

* 

22 

i 

23 

+ 

24 

.9923769 

25 

* 

2b 

STO  11 

27 

END 
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1 7 .       Change 

a.  PURPOSE 

This  subrouttne  ts  used  to  expedtte  the  changing  of  up 

to  fTve  of  the  Input  parameters  whenever  a  design  restraint 

ts  exceeded  and/or  at  the  end  of  the  main  program. 

b.  ASSUMPTIONS 

The  primary  parameters  which  will  require  changing  are: 

a.  Forward  velocity,  VF;     (Kts) 

b.  blade  twist,  TW;         (rads) 

c.  lift  curve  slope,  a;      (per  rad) 

d.  tip  velocity.  VT;        (ft/sec) 

e.  weight,  W;  (lbs) 

C.   EQUATIONS 

VF  (ft/sec)  -  VF  (Kts)  *  K68894 
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d.   FLOWCHART 


start 


clear 
flag  04 


yes  = 


no  =  0 


set  flags 
04.  2T 


yes 


no 


go  to  INO 
R31* 


stop 


DISPLAY: 


VF(Kts) 


TW 


VT 


prompt: 
VF(Kts) 


orompt: 
TWIST  -? 


change  Kts  to 
ft/sec  &  store 
a25 


prompt: 
vt  =  ? 
stnrfi  R 


prompt: 
a  -  ? 
store  R 


prompt: 
W  »  ? 


clear 
flag  27 

1 


NOTE:  ♦  returns  and 
reruns  mafn 
program  wfth 
new  data. 
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e.  PftOGHAMS    AND    SUflROUllNES    USED 

•    W8S    • 

f.  PftOGAAM   LISTING 


81*l8L  °CHG' 
62  °QGH° 
93  mSTO  31 
84  CF  94 

85*LBL  86 
96  "CHhHGE'^" 

97  PROHPT 

98  !^=9'^ 

99  GTO  97 
18  SF  94 
11  SF  27 

I?  lU   5  VT  H° 
13  PROHPT 


31*LBL  S 

33  PROHPT 

34  STO  13 

35  GTO  85 

36*LBL  E 

37  "U=^' 

38  PROHPT 

39  STO  19 

49*LBL  85 
a   CF  27 
42  GTO  86 


14*LBL  M 

16  PROHPT 

17  1.68894 

18  * 

1 Q  qj fi  OP 

29  GTO  95 


43*LBL  87 
44  FS'"  94 
45  GTO  IHB  ': 
46  END 


21*LBL  B 
12  •THIST='' 

23  PROMPT 

24  STO  29 

25  GTO  85 

26*lBL  C 

27     '3=":"' 

28  PROHPT 

29  STO  69 
39  GTO  85 
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D.        HP    41-CV    SAMPLE    OUTPUT 

Depanding   en      the   value     of    forward      valocity    •thrr<=-      ar 
threa   possible      *:ypes   of    outpur      that    tha    H?      41-C7   compu-a: 
program   is      capabls    of      producing.         Hxamplas      of   -ha      thrs 
oossibla    casas   are  listed    below. 


CASE    1    :      Forward   velocity,    V?,    is    l^ss    -han    s-all 
onset    velocitv. 


r5T_7q-:     7Ci  it;-; 
!^  i-=   7i  .  ■  -Ji  iJ--= 

PP=6. 986866 

VS= 137. 762848 
V6EK=7fe, 968756 
yi1fiK=  163. 443338 

CYCLIC=8. 666866 
CnLL= 17. 896967 

d278=4. 388796 

/9ei)=4. 388798 

PS=8. 966666 

PH=6. 868968 

PT=1..  816. 255686 
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CASE    2    :      Forward  velocity,    VF,    is   gr=ra 
stall    or.S€-    velocity. 


,__    -  V 


PI=113. 

37 

pn=352. 

43 

P?=794. 

24 

¥S=137. 

■b 

VBEK=76 

.86 

VHflX=lb 

3.45 

CYCLIC= 

-9.3? 

C0LL=2e 

.76 

i279=i7 

.75 

i90D=-l 

.93 

XS=8.74 

KS=i?,:*6 
PS=i48.98 

?«="'53  71 


PT=1 ,57"   16 
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CASZ    3    :      Forward  velocity,    7F,    is    greater    th; 
th9   maximum   forward   v^icci-.y    possib 


^0=368.92 
pp=844.72 

VBEK=7fc..S6 
V!1fiX=163.45 

VF  >  s'nm 

CtCLIC=-18.55 
0011=22.88 

i276=13.o7 

i98D=-i,23 

XS=8.73 

KS=i.86 

ps=i.j6.ie 

P?^=285.23 


FT=L76b.l5 
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APPENDIX    C 

A.        lEM    3  0  33    PROGRA.'I    DOCaMENT  ATI  ON 

This  program  calculates  the  power  required  -.  -  fly  a 
helicopter  in  forward,  straight  and  level  high  speed  fligh*. 
In  the  sections  tiat  follow  all  inpu-s  and  implement  a- ion 
requirements  are  specified.  The  problem  solving  Dne-.hodclcgy 
used    is    as    described    in    Chaorer    II. 
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B.        INPUT    DATA    REQUIRED 

Ir.  calculating  the  performance  of  a  aelicopi^r  it  is 
nacessary  to  defin?  a  set  cf  force  conditions,  environmental 
conditions,  and  physical  conditions.  The  force  conditions 
required  are  the  aircraft  gross  weigh-  (W)  in  pounds,  the 
.Tiaxinium  lift  coefficient  (CLMAX)  ,  and  tha  th-=  coefficient  of 
list  at  zero  angle  of  attack  (CDO) .  The  environmental 
conditons  include,  forward  velocity  (V?K|  in  Icncts,  -peed  of 
sound  (SV2L)  in  feet/sec,  rotor  height  above  the  ground  (H) 
in  feet,  and  air  density  (RHO)  in  Ib-ft  /sec  .  For  sake  of 
user  simplicity,  the  the  speed  of  scuni  and  air  density  ar= 
generated  within  the  prcgram  and  are  comparable  to  th=: 
values  found  in  a  standard  atmosphere  table.  Finally,  the 
physical  conditions  required  are  the  rotor  radius  (R)  in 
feet,  tip  velocity  (VT)  in  ft/sec,  numoer  of  blades  (b) , 
main  rotor  chord  (C)  in  feet,  flat  plate  area  (FPA)  in  ft  , 
geometric  twist  of  the  rotor  (TWIST)  in  radians,  and  airfoil 
lisft   curve   slcpe    (CLA)     in  per   radians. 
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C.        HELICOPTER    SAMPLE    DATA 

Tabls      C.I    below      illustrates  the      fDzmat  required      when 
inpating    data. 


TABLE    C.I        HELICOPTER     SAi^_?LE    IJIPaT    DATA 


RADIUS  CDO                                   W  VT  ?:PA 

22.  .01075                   10512.  738.0  17. 

TWIST  NO.     BLADES               CHORD  INITIAL    VFK         CLA 

-.17a5                     2.                        2-25  0.0  5.73 

H                             CLMAX  NO.    ENGINES  TYPE    AIRFOIL  DA 

1000.                   I.a                            2.  7.  1000 
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D.    IBM  PROGRAM  FLOWCHART 


Start 


read 

fnput. 

data 


calculate 
dfsk  area 


calculate 
solfdtty 


calculate 

induced 

velocity 


calculate 
AMAX 


calculate 

coef  of 

thrust 


calculate 
G.E,  ratto 


© 


calculate 

tTp-loss 

factor 


Calculate 
VMAX 


calculate 

fnduced 

power 


YES 


apply  C£ 
ratio 


GE 


1.0 


correct 

for  tfp 

losses 


calculate 

profile 

power 


3 


1  10 


YES 


® 


calculate 

parasite 

power 


calculate 

stall 
parameters 


calculate 
TUT2, 
T3,T4 


YES 


calculate 
A11.A12. 
A13.A14 


CPS  =  0 


calculate 
fnftfal 
lambda 


calculate 
XO,XS 


calculate 
Tnf low 
ratTo 


YES  ^z'^^.   NO 

l^VERR^ 


calculate 

cycHc, 
col lecttve 


A  =  A 


NO 

^-^ICav^lVs 

YES 

' 

^^"-^^^                1 

1 

calculate 
KS 

■ 

KS  -  1 

i 

1 

apply 

stall 

correctfon 

1 11 


calculate 
a270 


calculate 
comp  power 


calculate 
total 
power 


calculate 
VOPT 


set  VF 

VOPT 


NO 


'0=/ 


YES 


^ 

x270>AM? 

rv 

•i 

VF    =    VF    + 
7VF 

VF    » 
VF-VVF 

go  to  F 


stop 


6 


calculate 
a270 


6 
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2.        IBM    PHCGRAM    LISTING 

This      sectior.      contains      zh^    listing      of      the      IB'A      3033 
computsr    program   dsveloped  in   -his    report. 
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MAIN    ROTOR    HEIGHT    ABOVE    GROUND. 

AIR    DENSITY    (RhOi 

SONIC    VELOCITY ....• 

BLADE    GEOMETRIC    TWIST 

TYPE    AIRFCIL 

MAXIMUM  2-0  LIFT  COEFFICIENT. 
2-D  STATIC  STALL  ANGLE  (AMAX) 
LIFT   CLRVE    SLOPE    (/RAO) 

ZERO-LIFT  DRAG  CQEEFIC  I  ENT  .  .. 
CRITICAL  MACH  NC  (FCR  Ci  -  Oi 
INITIAL    FORWARD   VELOCITY    (XTi 


•  ,Gi2.o/, 
f  Fi2.o / r 


tG12.6/  f 
,Gi2.5// 
»Gi  2  .6/ f 
,Gi2.6/  , 
,G12.7/ , 
»Gi2.6  /  , 

I  F  i  2  •  6  /  f 

,Fi2.6// 

•       *    tGl  ^  aS  /  } 


END 
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?.        SAMPLE    OF    IBM   COMPUTEP   OaTPQT 

This  secticn  contains  an  example  run  of  -he  I3M  coDnpa-rr 
prograa,  utilizing  AH1-J  Cobra  data,  starting  a-  a  forward 
velocity    of    120    knots   and    terminatina   at    7:iAX. 
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hIGH    5P5EO    FORwARO    FLIGHT    ANALYSIS 
INPUT     DATA 

RADIUS • ZZ.QQOO 

MAIN    ROTOR    CHCRC...... • 2«25000 

NUMBER    OF    MAIN    PQTOR    BLADES. 2.00000 

AICRAFT    GROSS    wtlGHT.. 10612.0 

ROTOR    TIP    VELOCITY 738.000 

HORIZONTAL     FLAT    PLATE    AREA 17.0000 

NUMBER    OF    ENGINES    IN   HELICOPTER..  2.00000 

DENSITY    ALTITUDE 1000.00 

MAIN    ROTOR    hEIGhT    ABOVE    GROUND...  1000.00 

AIR    DENSITY    (RHCJ 0.002309 

SUNIC    VELOCITY 1113.04 

BLADE    GEOMETRIC    TWIST -.174500 

TYPE    AIRFOIL 7.00000 

MAXIMUM   2-D    LIFT    COEFFICIENT 1.40000 

2rQ, STATIC    STAL^    ANGLE     (AMAX) 2443280 

LIFT    CURVc    SlCPc    (/ftAO) 5.7j000 

ZERO-LIFT    CRAG    COEEFICIENT 0.010750 

CRITICAL    MACH    NC    IFOR    CL    =    0)..  0.720000 

INITIAL   FORUARO    VELOCITY    (KT)....       120. COO 

COEFFICIENT    GF    THRUST 554153E-02 

INDUCED    VELOCITY.. 33.8469 

DISC    AREA 1522. o2 

SOLIDI  TY 65  0197E-0i 

TIP-LOSS    FACTOR 947362 

GROUND    EFFECT    RATIO 1.00000 


FORWARD    VELOCITY     IN    KNOTS         =         120.  OCO 


PARASITE    DRAG    =       80t.258 

DISK    PLANE    ANGLE    CF    ATTACK    =    -.759760E-01 


DYNAMIC    PARAMETERS 

All 0.625400 

A  12 0.8069  2  6 

A13.... ....•••....•...•..  0.  5  733  3  8 

A 14..  ... 1.1754  3  5 

Tl 0.4676  0  2 

T2...... • 0.319142 

T3 0.2182  96 

T4. 0.12582  6 

INFLOW    kaTIC    i LAMBCAi  ...  . -0.040381 


1  2U 


STALL    POWER    C;iLCULAT  IONS 

RT -.lC57  78t-0i 

INBOARD    STALL    CQRRECTICN    FACTCR 0 

STALL    POWER    COEFFICIENT 0 


ANGLE    OF     ATTACH    CALCULATIONS 

LONGITUDINAL    CYCLIC   ANGLE.. -0.109613 

LONGITUDINAL    CCLLECTI Vc    ANGLE...  0.314954 

ALPHA(901      (JEG) -0.043036 

ALPHA<270J     (OEG) 12.511999 


HIGH    SPEcC    MACH    EFFECTS 

ADVANCING    ELADc    TIP   MACH    NUMBER...    .345135 

CRITICAL    ^ACH    NUMBER 7205^3 

DRAG    DIVERGENCE   MACH    NUMBER 6459i9£-01 

COMPRESSIBILITY    POUER    COEFF 52i491£-04 

PCWfcR   REQUIRED 

INDUCED    PCrtER    =  151.5^5 

PROFILE    PC^ER    =  29t),464 

PARASITE    POWEP    =  297.103 

COMPRESiaiLITV   POWER  =       134.001 

STALL   POWER    =  .0 

TOTAL    POWER    REQUIRED    =       879.112 

FCRWAPO    VELOCITY     IN    KNOTS         =         130.000 

^;c  4 4c 4C4C  4(  4c«  4:  4c  4  4  « 4( -< 4(4c-;i' 44c:te  «  «4:  4: 4i  «4c 4c  ««4: « -4e 4C4C ^ «  ^c 4(4c  4c  4<4c  4^ 4:  4u9c 4c 

PARASITE    DRAG    =       94^.234 

DISK    PLANE    ANGLE    OF    ATTACK    =    -.8916o4E-01 


OYNAI-IC    P^JRAMETERS 

All..... 0.6  8017  2 

A  12. 0.8803  7  7 

A  13 0.6  258  8  2 

Al>^...... 1.20  74  72 

Tl • 0.4  708  75 

T2 • ••..,••..  0.325344 

T3 0.221234 

T  4.. ....••. • 0.136  796 

INFLOW    RATIO    (LAMBDA)...... -0.044663 
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STALL    FC^ER    C^LCULATICNS 

Rl  •••••• .•...••••.••• -•  79341  c£>  0  2 

INBOARD    STALL    CORRECTICN    FACTCR 0 

STALL    POWEF    COEFFICIENT 0 


ANGLE    OF    ATTACK    CALCULATICNS 

LCMGITUOINAL    CYCLIC   ANGLE -0.122350 

LCNoITUDINAL    CCLLECTIVE    ANGLE...  0. 326189 

ALPHA(90)     (OEG).... ..•••  -0.291443 

AL?HAi270)     (DcCi. ..•••• 13.728103 


HIGH    SPEEC   MACH    EFFECTS 

ADVANCING    ELAOE    TIP   MACH    NUMBER...    .860209 
CRITICAL    MACh    ^UMBER 723294 

DRAG    DIVERGENCE   MACH    NUMBER.......    .770i57£-Oi 

COMPRESS! eicIT>    POWER    COEFF. ...... .    .630606E-Q4 

FCWER    REQUIRED 

3SC  4:  ^t^c^c^c  ^c^^c  4: 4<  44i«  «  44e4c  «  «4c4c  4(4:4:  4c4c«4(  4:4:«««  4e4(4c  4e4c  «4(4:  «  ^^^c  ««  44(4tJS'«^'.^ 

INDUCED    PCi^ER    =  139.900 

PROFILE    PC«£R    =  308.957 

PARASITE    PCWEF   =  377.740 

COMPRESiaiLITY   POWER  =      162.039 

STALL    POWER    =  .0 

TOTAL    POWER    RiQUi^EO    -      988. 63o 

4t  «***  *4:*  ***««««** 4=  4" *♦«**********»  4t4t**4'*«4c 4  4***44=* 

FCRi^ARO    VELOCITY    IM    KNOTS         =         140.000 

4c4c4[4::(c«4u«i:«c:4i:»^:$i;ji  4«:4e  4c  4**  ************  ^gc^c^ic  4c*  ****  44:**  4*  * 

PARASITE    DRAG    =       10^7.41 

DISK    PLANE    ANGLE    CF    ATTACK    =    -.103412 


OYNAyJC    P;aRAM£TERS 

All 0.7356  3  0 

A  12 • 0.  9  565  63 

A 13 • ••  0.6  79&5  9 

Al'f 1.2426  2  9 

Tl 0.4  74411 

T2 ••........ 0.3320  42 

T3 0.22440  7 

T4 0.1478  8  8 

INFLOW    RATIC     ( LAM8CA) -0.350098 
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STALS-    PG'^ER    C -^LC  ULATI  C.NS 
RT -.3a97  0I£-0i 

INBOARD    STALL    CQRRECTICN    FACTOR 0 

STALL     POVitfi    COEFFICIENT 0 


ANGLE    OF    ATTACK   CALCULATIONS 

LONGITUDINAL    CYCLIC    ANGLE. ♦...•.    -0.136576 
LQNGITUOIN^lL    CCLLECTIVE    ANGLE...       0.339933 

ALPHA(90}     (CEG).. -0.520492 

ALPHAI270)     (DEC) 15.129103 


HIGH    SPEEC   MACH    EFFECTS 

AOVANCIlMG    tLAOc    TIP   iMACH    NUMBER...    .8  7548  4 

CRITICAL    MACh    ^UMB£R... 725882 

DRAG    DIVERGENCE   MACH    NUMBER 896012E-01 

CQMPRESSieiLlTV    POk^ER    COEFF 745a73E-04 

PCWER    REQUIRED 

INDUCED     PC^vER    =  129.913 

PROFILE    PC^ER    =  322.449 

PARASITE    PCwEf    =  471.789 

COMPR£SIBILIT>    POKER  =      191.657 

STALL    POWER    =  .0 

TOTAL    POWER    REQUIRED    =       1115.81 

;;::»4c4c4c4c4c4c«4  444c«  44:4:4c  4>4c4c  4c4c;ic  4:4c  4c« 4c 4c4c4c4' 4c4c4c  4c4(4c4c4c  4c  4^i: 4c 4^  «-i;i4: 

FCRWARC    VELOCITY     IN    KNOTS         =         150.000 

4c4:jc4c:ic;:c;i^:Cc4c4i4c44c44c4c4c4>4:4c4c4c4:4c4c4c4:4c4c4c4c4c4c4c4c4c4c4>4c4c4c44c4c4c44:4c 

PARASI TE    DRAG    =       1259.78 

DISK    PLANE    ANGLE    OF    ATTACK    =    -.118713 


OYNAyiC    PARAMETERS 

All... • •••• 0.  791849 

A  12 .•• 1.  Q3'tl7C 

A  13 .  0.7348  01 

A14«. .•.......«•• • 1.281051 

Tl 0.4  782  0  8 

T2. 0.3  3923  7 

T3 0.227315 

T4 0.159103 

INFLOV^    RATIC    (  LAMBCAi -0.05o720 
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STALL    POWER    CALCULATIONS 
i^T 2307  7  7E-0  2 

VALUES    FOP    DETERMINING    STALL    COEFFICIENT 

GAMMA      =  -C. 26459350 

CS  =  -C.14755023 

BS  =  a. 3 244 962 7 

XS  =         0.7921403 

XO  =  1.C674362 

INBOARD    STALL    CORRECTION    FACTOR 66221dC94 

STALL    POWER    COEFFICIENT .  .     .3120259  SE-O'^ 


ANGLE    OF     ATTACK    CALCULATIONS 

LONGITUDINAL    CYCLIC   ANGLE -0.152545 

LCNGITUDINAL    CCLLECTIVE    ANGLE...  0.356377 

ALPHA{90)     (CEG  !•••  •••  •• -0.738637 

ALPHA{270i     (DECi.....  ••••••.....  16.740753 


HIGH    SPEEi    MACH    EFFECTS 

ADVANCING    ELADt   TIP   MACH    NUMBER...    .890o5a 

CRITICAL    MACH    NUMBER .728348 

DRAG    DIVERGENCE    MACH    NUMBER.. 102309 

COMPRESSIBILITY    POWER    CJEFF 867333E-04 

PCWER    RECUIRED 

INDUCED    PCy»ER    =  121.281 

PROFILE     PCiiER    =  336.941 

PARASITE    POWER   =  580.279 

COMPRESIBILITY    POWER  =      223.009 

STALL    POWER    =  80.1774 

TOTAL    POWES    REQUIRED    =      1341.69 

4c  4c 4c 4:4c  ^4c««  4c  4 44c 4c  44c4c4c4>4c4c  4c4c4c  4c4c4c4c4c4c4e4c4c  «4c4c  4c4c  44e4c  4>  44c4c4c  44c4c 

FCRWARO    VELOCITY     IN    KNOTS  =  160.000 

4c4c4c4c4c4e4c4c4t4  4  4c4c4c  444c  4  44c4c  «44  4c4i4c4c4c44c4  4c  4:-44c  4  4  444:  4c  44  44  «4c^ 

PARASITE    DRAG    =       1433.35 

DISK    PLANE    ANGLE    OF    ATTACK    =    -.135068 
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DYNAMIC    PARAMETERS 

All ••••••..••. 0.  848907 

A  12 i«  1133  9  7 

A  13 Q.  7<9144S 

A 14  ••....•«•••  ..•• 1.3228  9  9 

Tl •••• 0.4822  67 

T2 0.34692  7 

Ti 0.2  314  5  7 

T4..... 0.170433 

INFLOW    RATIO    i LAMfiOAi -0.064391 


STALL    POWER    CALCULATIONS 
RT...  ... 1 133 3 2E -01 

VALUES    FOP   OETERMIiNING    STALL    COEFFICIENT 

GAMMA      =  -C. 30194968 

CS  =  -C. 17313443 

BS  =  0.3  6384336 

XS  =         0.739  8  598 

XO  =         1*3366751 

INBOARD    STALL    CORRcCTICN    FACTOR. ...♦•       i.OCOOOCOO 
STALL    POWER    COEFFICIENT.... .  603492^  Ci-04 


ANGLE    Of    ATTACK    CALCULATIONS 

LONGITUDINAL    CYCLIC   ANGLE.. -0.17C333 

LONGITUDINAL    CCLLECTIVE    ANGLE...  0.375745 

ALPHA(90}      (OEGi.. -0.949361 

ALPHA(270)     (CEO 18.391431 


HIGH    SPEEC   MACH    EFFECTS 

ADVANCING    BLADE    TIP   M^CH    NUMBER...    .905332 

CRITICAL    MACh    NUMBER... .73C729 

DRAG    DIVERGENCE   MACH    NUMBER. .115103 

COMPRESSIBILITY    POWER    COEFF 997223E-04 

PCWER    RE  GUI  RED 

**:(i;^4c:ic:$c*4c***3;c4c:^  44::i^  ;;ic  4**  ****«**********<(*********  ******* 

INDUCED    PC'AER    =  113.722 

PROFILE     PC»«ER    =  352.433 

PARASITE    FGWEP    =  7C4.243 

COMPRESieiLlTV    POWER  =       25o.243 

STALL    POkEH    -  135. 38o 

TOTAL    PQwER    REQUIRED    =       1582.23 
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FCRwARJ    VELOCITY     IN    KNOTS         =         163.359 

***4c*  4c *«:**«*♦♦  44c4c4c44c*«4«4c4c4'4c4c4c4t4c4c4c4c*4c4t*4'*4c«44c«4:4'4c4t 

PARASITE    DRAG    =      1494.16 

DISK    PLANE    ANGLE    uF    ATTACK    =    -.1407S9 


OVNAMC    PARAMETERS 

All • 0.8682  7  7 

A  12 ••••• 1.  141191 

A 13 •••.•••••••...•...•.•...•  0. 810842 

Al-f • 1.  3377  5c 

Tl 0.4  836  8  9 

T2....». 0.349o22 

T3 0.2327  34 

T4 J.  1742  9  7 

INFLOW    RATIC    (LAHSCA) -0.067526 


STALL    PCV^ER    CALCULATICNS 
RT • 1559  8  7E-01 

VALUES    FOF    DETERMINING   STALL    COEFFICIENT 

GAMMA      =  -0.31570  089 

CS  =  -C. 18555188 

3S  =  0.36093823 

XS  -         0.7336488 

XO  =  1.^493771 

INBOARD    STALL    CORRECTION    FACTOR 1.00000000 

STALL     POWEF    COEFFICIENT.. .61  lCd76CE-04 


ANGLE    OF    ATTACK    CALCULATICNS 

LCNGITUOINAL    CYCLIC    ANGLE •  -0.177082 

LONGITUDINAL    CCLLECTIVE    ANGLE...  0.382947 

ALPHAi90i     (CEG) •«.. -1.019017 

ALPHA(270)     (OEG) 19.272018 


HIGH    SPEEC   MACH    EFFECTS 

ADVANCING    BLADE   TIP    MACH    NUMBER...    .91092  9 

CRITICAL    .MACh    ^UMBER 731516 

DRAG    DIVERGENCE   MACH    NUMBER. 119412 

COMPRESSieiLlTl    POrtER    COEFF. i04241£-0 
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PCWER    REQUIRED 

4c;«c4i4c:4cji::ic:ic;4c:»:«[4«:4i3ji44c4c34c«:p:9c  :«<  «V:;c4c4:jOc:tc««;Cc4c  4:4:^  ¥«  "(c^ic^c:}:  :$utc:ts:«:  :«:;^#:ic:(:^  4: 

INDUCED    PC*ER    =  111.338 

PROFILE    PC^ER    =  357.861 

PARASITE    FCWEP   ^  749.536 

COMPRESIBILITY    POWER  =       267.855 

STALL    POWER    =  157.023 

TOTAL    POWER    REQUIRED    =       1643.66 

INITIAL    STALL    CNSET    VELOCITY    APPROXIMATION    (KTS)     ^      76. 

DIFFERENCE    a£Tl.EEN    A270    AND   A.VAX    =       .2S9871E-02 
ACCEPATA8LE    ERi^OR    =       .lOOOOOE-02 

3(c:(c4c«:;c;^4c4c:ic«4«««  44c!«c«  ^^Cc^tc  :4i3t(;)e;ic«^4c«4c4c4c4c  3tc:(c4i4c  «  4««  :«c  :4i4c  «:4c  ;«l 

ANGLE    AT    STALL    ONSET=  13.9811 

STALL    GNSET    VELOCITY    =  131.913 

VELOCITY     HAX    E^DURANC£    =  76.0347 

MAXIMUM    FCPwARC   VELOCITY  -      163.359 

«  :^  4^4c 4c 4c4c^  «  :»  4 ««  4«4c  «  «4:«  «  «^  «$««  Xc  ««# 4c  «4c«  ««  «4::ic  Jtc  43(c  41:^ -^ 
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G.        COMPARISON    OF    PROGRAM    OUTPUT    VS    TEST  D AI A 

This    secticn    compares    the    output      of  the    HPU1-C7    and    13.1 

3033    computsr    programs      tc  actual   flight  test      data    gath^-re^ 
at    the    Naval    Weapons    System  Center. 


132 


POWER  REQUIRED  VS  VELOCITY 
FOR  THE  AH1-J  USING  THE  HP  41-CV 


TOTAL 


0.0  20.0  80.0  70. 0         100.0        123.0        100.0        176.0 

FORWARD  VELOCITY  (KT) 


Figure   c.1.      Powar    Curves  Generated   by    H?<41-CV   Prograi 
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POWER  REQUIRED  VS  VELOCITY 
FOR  THE  AH1-J  USING  THE  IB^  3033 
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Figure   c.  2.      Power  Curves   Generated  by   tlie   IBM   3033   Program. 
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POWER  HEOUIRED  ¥S  VELOCITY 


a  2 

U 


LEQEND 
TEST  FLSGHT  RESULTS 


O    IBM  3033  OUTPUT 
A    HP41-CV  OUTPUT 

___A_CCEPTA_B_LE_  ACCURAQY      TOTAL  POWER 

REQUIRED 


:x°y 


0.0  28. 0  SO.O  75.9  100.0  125.0  150.0 

FORWARD  VELOCITY  (XT) 


Figure  C- 3. 


Comparison    Between    Ccmputsr    Data   and    Actual 
Tesx  Plight    Data. 
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